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This  report  summarizes  the  results  of  a three-year  project  spon- 
sored by  the  U.S.  Army  Research  Office  and  directed  by  Dr.  James  A.  Smith, 
Principal  Investigator,  Colorado  State  University.  Dr.  D.  S.  Rimes 
received  the  Doctor  of  Philosophy  degree,  in  part,  for  the  work  reported 
here.  Mr.  K.  J.  Ranson,  Research  Associate,  Colorado  State  University, 
performed  many  of  the  analyses  and  coordinated  the  field  effort.  The 
major  objective  of  the  project  was  initially  to  apply  the  optical 
reflectance  modeling  procedures  (the  SRVC  model)  developed  previously 
by  the  authors  to  a forest  canopy  scene.  Subsequently,  project  objec- 
tives were  expanded  to  include  development  of  a thermal  canopy  exitance 
model.  This  latter  objective  was  made  possible  with  the  assistance  of 
personnel,  principally  Dr.  E.  Link,  of  the  U.S.  Army  Engineers  Waterways 
Experiment  Station.  An  extensive  data  base  was  obtained  in  cooperation 
with  WES  for  a lodgepole  pine  (Pinus  contort  a)  canopy  at  Leadville, 
Colorado  which  was  used  both  for  model  development  and  validation. 

Before  the  thermal  and  optical  terrain  feature  canopy  models  could  be 
applied,  techniques  for  describing  and  measuring  the  geometric  structure 
of  forest  canopies  had  to  be  ueveloped.  Rapid  optical  diffraction 
analysis  techniques  for  analyzing  ground  photographs  of  the  modeling 
trees  were  found  to  be  suitable.  The  details  of  the  models  developed 
and  their  applications  are  described  in  the  followin’  chapters  which 
contain  separate  introductions  and  conclusions.  Thi;  report  is  organ- 
ized as  follows: 


Chapter  2 describes  the  application  of  a modified  version  of  the 
solar  radiation  vegetation  canopy  (SRVC)  model  coupled  with  a numerical 
approach  to  estimate  solar  absorption  within  the  lodgepole  pine  canopy 
system.  The  Fortran  program  ABSORPT  performs  this  analysis  and  is 
presented  in  Appendix  B along  with  the  results  of  some  of  the  analyses 
mentioned  in  Chapter  2. 

Chapter  3 describes  the  TCSM  which  incorporates  the  information 
and  mathematics  derived  in  Chapter  3 along  with  the  thermal  radiant, 
convectional , and  t ranspirational  energy  exchanges.  The  model  predicts 
the  average  canopy  element  temperature  for  three  horizontal  canopy 
layers,  and  the  effective  radiant  temperature  above  and  within  the 
canopy  system.  The  results  of  some  of  the  analyses  mentioned  in  this 
chapter  and  the  Fortran  code  of  the  TCSM  are  presented  in  Appendix  C. 

Chapter  4 presents  the  summary  and  recommendations.  Appendix  A 
is  a reprint  of  the  abstracts  of  six  papers  or  reports  prepared  under 
partial  or  full  support  of  the  present  project. 
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Abstract 

The  absorption  of  total  and  spectral  solar  radiation  within 
vegetation  canopies  as  a function  of  solar  zenith  angle  needs  to  be 
quantitatively  described  for  agricultural,  ecological,  forestry,  and 
military  applications.  The  solar  radiation  vegetation  canopy  (SRVC) 
absorption  model  was  developed  to  physically  account  for  the  optical 
properties,  and  geometric  and  spatial  characteristics  of  canopy  ele- 
ments, and  direct  and  diffuse  components  of  irradiance.  Multiple  and 
directional  radiation  scattering  are  included.  The  model  predicts  the 
proportion  of  spectral  solar  absorption  in  three  horizontal  layers  and 
the  apparent  directional  reflectance  above  the  canopy.  Field  data 
were  collected  from  a cluster  of  four  lodgepole  pine  (Pinus  contor ta) 
trees.  Vertical  spectral  reflectance  above  the  canopy,  spectral  trans- 
mittance to  the  ground  layer,  geometric  measurements  of  canopy  ele- 
ments, and  optical  properties  of  canopy  elements  were  measured.  The 
model  was  then  applied  to  the  canopy,  and  the  reflectance  and  trans- 
mittance simulated  results  for  a theoretical  clear  day  were  compared 
with  the  measured  results.  The  simulated  results  showed  that  relative- 
ly large  differentials  occurred  in  spectral  absorption  by  canopy 
layers,  especially  in  the  photosynthetic  active  radiation  region  as  a 
function  of  solar  zenith  angle.  In  addition,  the  proportion  of  total 
global  irradiance  absorbed  by  individual  layers  varied  as  a function 
of  solar  zenith  angle.  However,  the  proportion  of  both  total  and 


spectral  global  irradiance  absorbed  by  the  entire  canopy  system  was 
relatively  constant  with  solar  zenith  angle. 

Introduction 

The  manner  in  which  a vegetation  canopy  absorbs  solar  radiation 
has  an  important  effect  on  the  thermal  properties  of  the  canopy  and 
the  photosynthetic  efficiency  of  the  canopy.  Thus,  an  understanding 
of  these  principles  is  important  in  remote  sensing  with  respect  to 
military,  agricultural,  forestry,  and  ecological  applications.  For 
example,  in  recent  years  the  thermal  region  of  the  electromagnetic 
spectrum  has  received  keen  interest  in  the  remote  sensing  field.  This 
region  may  add  valuable  additional  information  to  make  inferences 
concerning  the  characteristics  of  vegetation  canopies.  However, 
before  the  thermal  emission  characteristics  of  a canopy  can  be  under- 
stood, the  manner  in  which  the  canopy  absorbs  solar  radiation  must  be 
studied.  In  the  field  of  agriculture  there  is  strong  evidence  that 
sclar  radiation  distribution  within  a canopy  as  a function  of  canopy 
structure  strongly  affects  the  productivity  of  the  canopy  (Vidovic, 

1973;  Rhodes,  1971;  and  Donald,  1961). 

Physically  based  mathematical  models  serve  as  convenient  tools 
in  studying  the  complex  radiation-vegetation  interactions.  The 
objective  of  the  study  was  to  develop  a mathematical  physically  based 
model  to  study  the  manner  in  which  spectral  and  total  solar  radiation 
as  a function  of  solar  zenith  angle  are  absorbed  in  vegetation  canopies. 
The  following  describes  the  absorption  model  and  the  application  of 
the  model  to  a lodgepole  pine  (P i mis  contorta)  canopy  at  Lcadville, 
Colorado  for  which  a unique  data  base  was  collected  during  1977.  A 


complete  study  site  description  is  given  by  Ranson,  Kirchner,  and 

Smith  (1978).  The  specific  canopy  modeled  consisted  of  a cluster  of 

4 lodgepole  pine  trees  with  the  mean  statistics:  6.0  m height,  30 

yr.  age,  13.2  cm  diameter  breast  height,  and  a surrounding  stand  of 
2 , 

102  m /hectare  basal  area. 

Solar  Radiation  Canopy  Models 

Several  deterministic  models  have  been  developed  to  study  the 
interactions  of  solar  radiation  within  vegetation  canopies.  Allen 
and  Richardson  (1968),  Alderfer  and  Gates  (1971),  and  Suits  (1972) 
have  adapted  a system  of  simultaneous  differential  equations,  developed 
by  Kubelka  and  Munk  (1931),  in  various  ways  to  vegetation  canopies. 
Suits  (1972)  developed  a model  which  includes  geometric  effects  and 
predicts  non-Lambert ian  characteristics  of  vegetation  canopies. 

Chance  and  LeMaster  (1978)  have  derived  a light  absorption  model  for 
vegetative  plant  canopies  from  the  Suits  reflectance  model  (1972). 

Another  approach  developed  by  Oliver  and  Smith  (1974)  is  the 
solar  radiation  vegetation  canopy  (SRVC)  model.  This  model  simulates 
the  solar  radiation  flow  through  the  canopy  by  utilizing  physical 
laws  and  Monte  Carlo  techniques.  This  stochastic  model  originally 
predicted  the  diurnal  apparent  directional  spectral  reflectance  of 
a vegetation  canopy. 

It  is  believed  that  the  ray  tracing  technique  utilized  in  the 
SRVC  approach  is  advantageous  as  applied  to  solar  radiation  inter- 
actions within  vegetation  canopies  for  several  reasons.  The  total 
effect  of  all  possible  events  can  be  simulated  if  one  knows  the 
probability  for  each  step  in  a sequence  of  events.  Thus,  as  new 
knowledge  becomes  available  on  the  probabilities  for  each  step,  the 
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SRVC  framework  can  readily  accept  it.  For  example,  if  a researcher 
describes  a non-Lambert ian  reflection  or  transmission  distribution  as 
a function  of  the  source  direction  and  leaf  orientation,  this  infor- 
mation could  be  incorporated  relatively  easily  within  the  SRVC  frame- 
work. This  framework  has  other  advantages  when  applied  to  vegetation 
canopy  modeling: 

(1)  Such  a general  framework  can  be  easily  modified  to  include  addi- 
tional considerations  without  having  to  examine  their  effect  on 
the  solution  to  differential  equations  as  in  the  deterministic 
models . 

(2)  The  model  can  be  modified  to  accept  any  reasonable  number  of 
components  within  a scene.  Thus,  one  could  model  a scene  as 
complex  as  time  permits  to  obtain  reasonable  geometric  and 
spectral  data  of  the  components. 

(3)  Any  relevant  parameter,  such  as  number  of  components,  type  of 
component,  their  reflectance,  transmittance,  and  absorptance 
angle  distributions,  surface  area,  and  spatial  dispersion,  may 
be  varied  in  any  desirable  fashion,  and  the  model  can  accept 
this  information. 

(4)  The  model  can  be  modified  to  any  reasonable  number  of  discretized 
inclination  and  azimuthal  angles  for  simulation. 

(5)  Diffuse  skylight  is  treated  as  a set  of  independent  source 
vectors. 

(6)  The  model  accounts  for  multiple  reflection  and  transmission  in 
both  upward  and  downward  directions. 


The  disadvantage  of  the  proposed  model  is  mainly  one  of  the  relatively 
large  computer  time  involved  per  run. 
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For  these  reasons,  the  SRVC  model  was  modified  to  produce  a 
version  of  the  model  to  study  the  solar  absorption  within  vegetation 
canopies.  A complete  description  of  the  original  SRVC  model  is 
presented  by  Oliver  and  Smith  (1974). 

SRVC  Absorption  Model 

An  abbreviated  las  at  ip t ion  of  the  SRVC  absorption  model  is  as 
follows.  The  SRVC  absorption  model  assumes  that  a vegetation  canopy 
is  composed  of  non-homogeneous  layers  of  Lambertian  elements  of  known 
geometric  arrangement,  statistical  composition,  and  optical  proper- 
ties. The  global  radiation  is  composed  of  direct  and  diffuse  sky 
radiation.  The  direct  solar  radiation  is  treated  as  a point  source, 
and  the  diffuse  radiation  is  divided  into  source  sectors.  These  source 
sectors  are  created  by  dividing  the  hemisphere  into  inclination  bands 
and  then  further  dividing  each  band  into  sectors  (Figure  1).  This 
spherical  coordinate  framework  serves  as  an  accounting  method  for 
radiation  transfer  above  and  within  the  canopy  system.  The  flux  from 
each  sector  is  simulated  as  a source  vector  (Figure  2).  The  inter- 
action of  each  initial  source  vector  from  the  sky  with  the  canopy  is 
then  calculated  independently.  The  model  utilizes  probabilities  which 
govern  the  distribution  of  gaps  within  the  vegetation  to  determine 
the  transition  of  source  vectors  from  point  to  point  within  the  canopy. 

The  formulation  developed  by  Idso  and  deWit  (1970),  which  is  a 
function  of  the  canopy's  geometry,  has  been  incorporated  to  predict 
the  probability  of  gap  in  the  direction  of  the  nine  hemispher ical 
bands  for  each  canopy  layer.  In  this  particular  study,  three  canopy 
layers  of  equal  height  were  defined  (Figure  2).  The  positive  binomial 
distribution  is  used  to  describe  these  probabilities.  Azimuthal 


r SECTION 

:CR  BAND  3 


DIRECTION 
f OR  BAND  I 


Horizontal  view 


NORTH 


SECTOR 


Figure  1 


Horizontal  and  vertical  views  of  the  9 hemispherical 
inclination  bands  which  are  divided  into  18  equal 
azimuthal  sectors.  If  one  rotates  the  horizontal  view 
about  the  Z axis  the  bands  would  occur  in  three-space. 
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Figure  2.  Schematic  of  a plant  canopy  approximated  by  stratified 
vegetation  layers  containing  a statistical  ensemble  of 
Lambertian  surfaces.  Two  interactions  are  shown  with 
the  corresponding  source  vectors  representing  reflection 
and  transmission  in  a Lambertian  manner. 
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symmetry  is  assumed.  The  probability  of  gap  is  equal  to  the  ratio  of 
the  projection  of  canopy  elements  in  any  particular  layer  to  the 
projection  of  the  underlying  soil  surface  for  each  hemispherical  band. 
For  a hemispherical  band  direction  0 (inclination  angle)  the  equation 
is : 

T LAI 

PGAP (?)  = 1 - S-g(O)  s'  ; PHIT (0)  = 1 - PGAP (0)  (1) 

sin (0) 

where : 

PGAP(0)  = probability  of  gap  in  direction  of  hemispherical 
band  0 

PHIT(O)  = probability  of  hit  in  direction  of  hemispherical 
band  0 

g(0)  = mean  canopy  projection  in  the  direction  of  hemi- 

spherical band  0 

LAI  = leaf  area  index 

S = index  of  spatial  dispersion. 

The  function  g(0)  is  determined  from  element  inclination  angle  distri- 
butions which  describe  the  orientation  of  the  elements  in  a canopy 
layer.  The  derivation  and  computational  procedure  is  presented  by 
deWit  (1965).  The  parameter  S ranges  from  0 to  1 and  is  an  index  of 
denseness  or  spatial  dispersion  of  the  components  in  a canopy  layer. 

As  S approaches  1,  the  more  regular  the  dispersion  of  components  and 
the  less  frequent  a gap  is  encountered.  The  leaf  area  index  (LAI)  of  a 
canopy  layer  is  equal  to  the  ratio  of  the  total  one-sided  element 
area  within  a layer  to  the  area  of  the  underlying  soil  area.  For  a 
more  in-depth  discussion  of  the  above  theory  and  the  required 
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measurements  see  Kiraes,  Smith,  and  Berry  (1°78);  Smith  and  Berry 
(1976);  and  Idso  and  deWit  (1970). 

Wien  a canopy  element  or  a ground  element  is  hit,  a proportion 
of  the  flux  vector  is  reflected,  transmitted,  and  absorbed  into  a 
number  of  flux  vectors  which  simulate  a Lambertian  response  (Figure  2). 
The  direction  of  these  vectors  is  determined  by  the  element's  orien- 
tation. The  proportion  of  reflected  and  transmitted  flux  is  deter- 
mined by  the  spectral  characteristics  of  the  canopy  elements.  These 
resulting  flux  vectors  are  further  processed  in  a similar  fashion 
until  all  vectors  are  essentially  zero,  indicating  absorption  by 
canopy  elements  and  ground,  or  escape  from  the  canopy. 

The  SRVC  model  predicts  the  apparent  directional  spectral  reflec- 
tance of  a canopy  in  the  nine  hemispherical  inclination  bands. 

The  SKVC  model  was  modified  in  the  following  manner  to  predict 
solar  absorption.  For  each  source  vector-element  interaction  we  know 
that : 


PF.  = (aX  + PX  + TX  )‘PE 
X X 


wnere : 


= the  proportion  of  spectral  solar  irradiance  (E  ) of  a 

e» 


wavelength  band  represented  by  a source  vector  incident 
on  a canopy  element 

a ^ = spectral  absorption  coefficient  of  canopy  clement 
0^  = spectral  reflection  coefficient  of  canopy  element 

= spectral  transmission  coefficient  of  canopy  element. 

The  reflectance  and  transmission  values  of  the  canopy  elements 
for  each  discrete  wavelength  are  input  to  the  model.  These  coeffi- 
cients are  different  for  different  types  of  canopy  elements,  and  thus 
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a. , p.,  x are  dependent  on  the  material  type  hit.  The  proportion  of 
AAA 

solar  irradiance  in  a discrete  wavelength  that  is  absorbed  by  each 
canopy  layer  at  a particular  solar  zenith  angle  is  obtained  by  s.  ::r.ing 
the  absorbed  proportion  of  all  the  source  vectors  incident  on  the 
canopy  eler.ients  occurring  within  a particular  layer. 

where : 

ct.  . = the  proportion  of  spectral  solar  irradiance  at 

wavelength  X within  layer  i,  i = 1,  2,  3,  4 where 
layer  4 designates  the  ground  at  solar  zenith  angle  z. 
The  SRVC  absorption  model  conserves  energy,  i.e.,  the  sum  of  canopy 
reflectance  and  absorptar.ee  equals  the  incoming  solar  irradiance. 

To  estimate  the  absorbed  total  solar  irradiance,  it  is  necessary 
to  numerically  integrate  the  spectral  absorption  at  a number  of 
discrete  wavelength  bands  over  the  entire  electromagnetic  spectrum. 
Because  computer  time  required  is  directly  proportional  to  the  number 
of  wavelengths  simulated  and  the  number  of  canopy  elements,  a simplified 
numerical  integration  scheme  was  devised  as  explained  later. 

Approach 

To  estimate  the  total  and  spectral  global  solar  irradiance  absorbed 
by  the  lodgepole  pine  canopy,  the  following  information  is  required: 
canopy  element  area  index  for  each  camoy  layer;  normalized  spectral 
solar  irradiance  curves  ranging  from  0°  - 85°  solar  zenith  angles; 
proportion  of  spectral  irradiance  absorbed  by  each  canopy  layer;  and 
the  mean  canopy  element  spectral  reflectance  and  transmittance  for 
each  discretized  wavelength.  The  entire  numerical  approach  is 


. i,z 


I 


all  interactions 
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summarized  in  Figure  3.  First  the  numerical  approach  will  be  presented 
followed  by  the  data  acquisition. 

Using  the  above  inforraation  and  interpolation  and  integration 
techniques,  the  absorbed  total  solar  radiation  within  each  canopy 
layer  is  estimated  as: 


where : 


I.  = approximated  integral  of  absorbed  total  solar  flux 
1 , z 

within  layer  i at  solar  zenith  angle  z 

= total  global  solar  irradiance  at  zenith  angle  z 

f,  = normalized  spectral  solar  irradiance  curve  at  zenith 
X , z 1 


angle  z 

a . = simulated  proportion  of  spectral  global  solar  irradiance 

A , 1 , Z 

absorbed  by  canopy  layer  i for  a mean  canopy  element 
reflectance  corresponding  to  wavelength  band  X and 
solar  zenith  angle  z. 

The  mean  absorbed  total  solar  flux  absorbed  per  unit  of  canopy 
element  surface  area  within  layer  i at  solar  zenith  angle  z can  then 
be  calculated  as: 


I. 

1 , Z _L  ,J! 

M2  = LA1i 

where: 

<!> 

i , z = average  absorbed  total  solar  flux  absorbed  in  layer  i 
M2 

at  solar  zenith  angle  z per  unit  area  of  element 
LAI^  = the  total  element  surface  area  within  layer  i. 


»• 
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The  spectral  global  irradiance  absorbed  by  the  canopy  is  estimated 
by  the  simulated  a,  . coefficients. 

A , 1 , z 

The  total  absorption  coefficient  of  a single  isolated  mean  canopy 
element  is  estimated  by: 


A 

Z 


■'VdA 


where : 


Az  = the  total  absorption  coefficient  of  a single  isolated  mean 
canopy  element  under  spectral  irradiance  conditions 
defined  by  solar  zenith  angle  z 
= spectral  absorptance  curve  as  a function  of  wavelength  of 
the  mean  canopy  element. 

The  method  of  numerical  integration  and  normalization  utilized 
throughout  this  study  is  as  follows.  The  problem  is  to  evaluate 

rb 

I = J f (x)dx 

when  f (x)  is  known  only  at  a finite  number  of  points.  First,  f (x)  is 
approximated  by  a polynomial  p(x)  and  then  p(x)  is  integrated  to 
obtain  I.  P(x)  is  described  by  using  Newton's  forward-difference 
interpolating  polynomial  (Conte  and  deBoor,  1965).  The  data  utilized 
had  non-uniform  discretization  intervals  and  a few  relatively  large 
intervals.  In  addition,  the  field  data  were  not  necessarily  smooth  by 
nature.  Under  these  conditions,  a higher  order  polynomial  fit  will 
not  necessarily  yield  a more  accurate  approximation  to  the  integral. 
Thus,  for  the  sake  of  simplicity  and  computing  time,  a first  degree 
polynomial  was  decided  to  be  adequate  for  the  purpose  of  integration. 
Thus,  the  composite  trapezoidal  rule  (Conte  and  deBoor,  1965)  was 
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employed  to  approximate  the  integral  of  the  above  curves  over  any 
desirable  [a,b]. 

The  n 'eessary  parameters  to  determine  the  total  and  global 
absorbed  radiation  were  determined  ns  follows.  The  total  canopy  element 
area  index  (LAI)  can  be  defined  as  the  total  surface  area  of  the  canopy 
elements  (e.g.,  leaves,  stems,  and  reproductive  structures)  divided  by 
the  projected  ground  area.  The  LAI  was  estimated  by  combining  the 
canopy  element  area  index  for  branches  (BAT)  and  needles  (NAI)  for  each 
canopy  layer.  The  procedure  involved  measuring  all  branch  diameters 
for  all  four  modeling  trees  near  the  bole  of  the  tree.  Regression  equa- 
tions for  lodgepole  pine  developed  by  Gary  (1976)  which  relate  branch 
diameter  to  the  total  l ranch  and  needle  surface  area  for  the  top,  middle, 
and  base  of  the  crown  were  utilised  to  derive  LAI,  BAI,  and  NAI  for  the 
total  tree.  The  proportions  of  NAI  and  BAI  for  each  canopy  layer  of 
equal  height  were  derived  from  Gary  (1976)  and  are  presented  in  Table  1. 
The  final  estimated  LAI  for  Layers  1,  2,  and  3 were  2.4,  4.5,  and  1.6, 
respectively,  for  a total  canopy  LAI  of  8.5. 

The  distribution  of  global  and  diffuse  solar  energy  as  a function 
of  wavelength  was  measured  on  a clear  day  at  the  study  site  using  the 
U.S.  Forest  Service  Circular  Variable  Filter  Spectrometer  (CVFS) . 

The  CVFS  uses  two  continuously  varying  interference  filters  for 
spectral  separation  in  the  visible  and  near  infrared  regions  with  a 
half  bandwidth  of  about  15-22  nm.  The  system  has  a 50  mm  camera  lens 
and  an  acrylic  diffuser  for  input  optics  and  a silicon  diode  operating 
in  the  photo  voltaic  mode  as  a detector.  A digital  readout  is  avail- 
able (letter  dated  14  December  1978  from  Robert  W.  Dana,  Physicist, 
Resources  Evaluation  Techniques  Program,  Rocky  Mountain  Forest  and 
Range  Experiment  Station,  U.S.  Forest  and  Range  Experiment  Station, 


i 


l 


% 


2-14 


Table  1.  Proportion  of  BAI  and  N'AI  of  respective  totals  for  the 
three  canopy  layers  of  equal  height.  Data  were  derived 
from  measurements  on  a single  lodgepole  pine  tree  (13.2  m 
height,  13  cm  DBH)  conducted  by  Gary  (1976). 


Proportion 


Layer  NAI  BAI 


1 

0.28 

0.23 

2 

0.53 

0.49 

3 

0.19 

0.28 

total 

1.00 

1.00 

I 
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U.S.  Forest  Service,  240  W.  Prospect  Street,  Fort  Collins,  Colorado 
80521).  Thirty-one  discrete  spectral  measurements  for  both  diffuse 
and  direct  radiation  were  recorded  within  the  range  of  .44  - 1.0  pm 
at  0S50,  1025,  and  1152  (Standard  Time)  during  August  6,  1976.  All 
measurements  were  taken  within  13  minutes  of  the  above  times.  These 
times  correspond  to  solar  zenith  angles  of  22°,  30°,  and  47°,  respec- 
tively. The  proportion  of  total  solar  irradiance  represented  in  the 
UV  range  (<.44  um)  and  two  IR  ranges  (1.0  - 2.0,  <2.0  pm)  were  esti- 
mated from  tables  presented  by  Kondrat'yev  (1965)  which  typify  a 
theoretical  clear  and  dry  atmosphere.  The  general  trends  of  the 
field  data  were  compared  with  those  theoretical  curves  (Kondrat'yev, 
1965)  by  normalizing  all  solar  spectral  curves  so  that  the  integral 
of  each  curve  equaled  1.0. 

This  computation  aided  in  the  ease  of  comparison  between  shifts 
in  the  above  spectral  curves  and  in  later  analysis.  The  computation 
involved  integrating  the  relative  magnitudes  of  the  points  within 
the  .44  - 1.0  pm  regions;  and  from  tabular  values  (Kondrat'yev,  1965), 
the  corresponding  proportion  of  total  solar  irradiance  in  the 
remaining  UV  and  IR  regions  were  obtained.  The  normalized  f^  values 
within  interval  (.44,  1.0)  were  then  calculated  as: 

, V<1-P> 

Ei  = ~ T7 

where : 

I 

f^  = normalized  f values 

Ei  = = 2 N)  whore  e (.44  - 1.0) 
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I*  = approximate  integral  on  interval  (.44,  1.0) 
p = proportion  of  total  solar  energy  outside  the  .44-1.0  ym 
region. 

The  measured  solar  spectrum  does  not  span  the  entire  solar  zenith 
angle  range.  However,  the  normalized  spectral  curves  as  measured  in 
the  field  were  very  similar  in  regard  to  spectral  trends  of  the 
appropriate  theoretical  spectral  curve  (Figure  4).  Further,  because 
of  the  relatively  small  change  in  atmospheric  path  length  between  the 
three  measured  curves  and  the  theoretical  curve,  the  trends  and  absolute 
normalized  magnitudes  are  similar.  In  addition,  these  curves  corres- 
pond relatively  well  to  the  spectral  irradiance  curves  reported  by 
Gates  (1966).  As  a consequence,  the  4 theoretical  curves  as  presented 
by  Kondrat'yev  (1965)  for  zenith  angles  of  0,  70,  80,  and  85°  (respec- 
tive atmospheric  path  lengths  of  1.0,  3.0,  6.0,  10.0)  were  utilized 
exclusively  in  this  study  for  consistency  (Figure  5).  The  proportion 
of  solar  irradiance  in  the  UV  and  IR  bands  for  the  four  theoretical 
curves  are  presented  in  Table  2.  These  four  curves  demonstrate  that 
for  a clear  and  dry  atmosphere,  as  zenith  angle  increases,  the  path 
length  through  the  atmosphere  and  the  shift  in  the  direction  of 
longer  wavelength  increases.  The  general  trends  in  spectral  shifts  of 
these  curves  were  believed  adequate. 

These  4 normalized  spectral  solar  curves  can  then  be  linearly 

interpolated  for  any  specific  zenith  angle  and  for  specified  discrete 

wavelength  intervals  to  obtain  the  desired  discretized  spectral  solar 

irradiance  function  (f.  ) for  zenith  angle  z and  at  wavelengths 

a , z 

A . , i = 1,34. 
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Figure  5.  Theoretical  solar  irradiance  curves  for  a clear  and  dry 
atmosphere  (from  Kondrat’yev,  1965). 





2-19 


Table  2.  Theoretical  proportion  of  total  solar  irradiance  In  the 
UV  and  IR  regions  as  a function  of  zenith  angle  (from 
Kondrat'yev,  1965). 


\;:m 

Solar 

zenith  angle 

0° 

70° 

O 

O 

85° 

<.44 

.089 

.048 

.021 

.007 

1.-2. 

.239 

.269 

.305 

.343 

>2 . 

.064 

.072 

.033 

.094 
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The  proportion  of  spectral  solar  irradiance  as  a function  of 
solar  zenith  angle  absorbed  by  each  canopy  layer  is  estimated  by  the 
SRVC  absorption  model.  To  simulate  all  permutations  for  each  of  the 
34  discretized  wavelengths,  2 canopy  components  (branches  and  needles), 
the  entire  range  of  solar  zenith  angles,  and  all  possible  diffuse/direct 
ratios  of  solar  irradiance  would  be  extremely  costly  and  thus  the 
following  approach  was  employed  to  reduce  the  number  of  required 
simulations. 

The  spectral  diffuse/direct  ratio  as  derived  from  the  U.S.  Forest 
Service  CVFR  data  collected  at  Leadville,  Colorado,  August  6,  1976  is 
presented  in  Figure  6.  The  ratio  increased  with  increasing  frequency 
and  increasing  solar  zenith  angle.  Wavelengths  near  the  UV  region  had 
a ratio  of  less  than  0.15.  A sensitivity  analysis  of  the  change  in 
absorption  coefficients  for  each  layer  versus  a change  in  the  diffuse/ 
direct  ratio  for  wavelengths  of  high  and  low  canopy  reflectance  and 
zenith  angles  of  0°  and  12°  was  completed.  In  the  ratio  range  of  0.0- 
0.2,  the  change  in  the  absorption  coefficient  is  less  than  0.03  in  all 
cases.  Thus,  it  is  believed  that  for  clear  sky  conditions  at 
Leadville,  Colorado,  a constant  diffuse/direct  ratio  of  0.06  is 
adequate  for  all  wavelengths. 

Rather  than  simulate  two  canopy  elements  (needles  and  branches) 
in  each  layer,  one  average  element  in  spectral  characteristics  was 
simulated  for  each  layer.  The  spectral  reflectance  of  the  branches 
and  needles  and  needle  transmittance  were  initially  measured  using  an 
Isco  Model  SR  Spectral  Radiometer.  Branches,  mosaics  of  needles,  and 
a barium  sulfate  reference  panel  were  utilized  to  obtain  reflectances. 

A metal  plate  with  a thin  slit  was  used  to  measure  needle  transmittance. 


Figure  6.  Measured  dif fuse/direct  ratios  of  solar  irradiance  as  a function  of  wavelength 
and  zenith  angle  (Z) . 
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In  addition,  a Modified  Barnes  Spectral  Master  Research  Radiometer  was 
utilized  to  measure  the  spectral  reflectance  curve  of  needles  and 
branches.  This  radiometer  had  a 2°  field  of  view  (no  needle  mosaics 
were  required)  and  produced  a continuous  curve  from  0.25  to  1.20  pm. 
Using  the  above  two  instruments,  14  needle  reflectance,  12  needle 
transmittance,  and  11  branch  reflectance  samples  were  taken  from 
various  portions  of  a tree.  The  respective  average  spectral  curves 
were  then  weighted  according  to  the  proportion  of  BAI  and  NAI  for 
each  canopy  layer  as  measured  by  Gary  (1976) (Table  1).  The  resulting 
spectral  reflectance  and  transmittance  curves  of  the  average  canopy 
element  for  each  layer  were  very  similar,  thus  only  one  set  of  mean 
spectral  reflectance  and  transmittance  curves  were  utilized  for  all 
three  canopy  layers  as  presented  in  Figure  7.  The  mean  spectral 
reflectance  curve  corresponded  closely  to  spectral  reflectance  curves 
of  Pinus  resinosa  needles  as  presented  by  F.gan  (1970)  . The  curves  in 
Figure  7 are  used  to  calculate  the  spectral  absorptance  curve  of  the 
average  canopy  element  for  all  three  layers  which  is  denoted  as  a . 

A 

It  was  also  essential  that  canopy  absorption  be  shown  insensitive 
to  small  changes  in  the  average  element  transmission.  The  proportions 
of  absorbed  solar  irradiance  for  a wavelength  of  high  canopy  reflec- 
tance (p  = .50  for  all  canopy  elements  including  the  ground)  for 
various  element  transmission  coefficients  and  for  zenith  angles  of 
0°  and  72°  were  simulated.  The  maximum  standard  deviation  of  the 
measured  needle  transmission  for  the  six  discrete  wavelengths  was  0.02 
and  corresponded  to  a near  infrared  wavelength  (Figure  7).  In  the 
above  sensitivity  analysis  a *.0.04  (two  standard  deviations)  trans- 
mission change  produced  a maximum  absorption  coefficient  change  on 
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the  order  of  0.01  for  Layer  1 and  a 0°  solar  zenith  angle.  The  stand- 
ard deviation  for  the  0.68  inn  band  was  0.0008.  At  low  reflecting 
wavelengths  any  reasonable  error  in  transmission  will  produce  a small 
change  in  the  absorption  coefficients.  In  addition,  it  was  found  that 

the  average  canopy  reflectance  was  correlated  with  the  average  canopy 
2 

transmission  (R“  = 0.62  with  reflectance  being  the  independent 
variable) . 

The  above  suggests  that  only  a fewT  mean  canopy  element  reflectance 
(MCR)  values  need  to  be  simulated  with  the  SRVC  absorption  model. 

Six  MCR's  were  simulated  ranging  from  0.0  to  0.S0  reflectance.  Each 
MCR  value  simulated  for  all  three  canopy  layers  has  a corresponding 
mean  canopy  element  transmission  factor  and  a ground  reflectance 
factor.  Thus,  the  MCR  defines  the  optical  properties  of  the  entire 
canopy  system.  The  above  insensitivity  and  correlation  of  transmission 
allows  one  to  make  the  transmission  factors  dependent  only  on  MCR 
rather  than  wavelength.  The  background,  composed  of  grasses  and 
litter,  was  similar  to  the  average  canopy  element  reflectance  except 
in  the  near  IR  region  where  the  background  was  consistently  lower  in 
rj  reflectance.  As  a consequence,  at  high  MCR  the  corresponding  reflec- 

tance of  the  background  was  scaled  down  relative  to  the  MCR.  The 
simulated  MCR  along  with  the  corresponding  canopy  element  transmission 
and  the  ground  reflectance  utilized  are  presented  in  Table  3. 

Four  runs  utilizing  the  SRVC  absorption  model  were  made  for  solar 
zenith  angles  of  0°,  47°,  72°,  and  87°  using  the  following  input 
parameters.  The  diffuse/direct  ratio  utilized  for  all  runs  was  0.06. 
The  canopy  clement  inclination  angle  distributions  were  derived  by 
the  laser  diffraction  technique  as  presented  by  Kimes,  Smith  and  Berry 

i 

i 


(1979).  The  S parameter  was  measured  to  be  0.1  by  measuring  the 
frequency  of  gap  using  black  and  white  photographs  which  were  taken 
looking  vertically  up  through  the  canopy.  A dot  grid  was  applied  to 
the  photographs  to  estimate  the  probability  of  gap.  The  S parameter 
was  then  derived  from  Equation  1.  The  LAI  values  for  each  layer  were 
derived  as  discussed  above. 

The  resulting  runs,  which  are  presented  in  the  Results  and 

Discussion  section,  can  be  linearly  interpolated  for  any  solar  zenith 

angle  and  MCR  which  corresponds  to  a specific  wavelength  (Figure  7) 

to  estimate  a.  . 

\,i,z 

The  total  global  solar  irradiance  at  zenith  angle  z was  measured 
using  a MARK  I-G  SOL-A-METER  Silicon  Cell  Pyranometer.  This  measure- 
ment, denoted  as  , when  multiplied  by  the  discretized  f^  ^ determines 
the  magnitude  of  the  spectral  curve,  whereas  f^  determines  the  shape 
of  the  curve.  Some  error  is  introduced  due  to  changing  solar  spectral 
effects  as  a function  of  solar  zenith  angle  and  the  sensitivity  of 
the  silicon  detector.  However,  this  will  effect  only  the  magnitude 
of  the  spectral  curve  and  not  the  theoretical  spectral  trend. 

Model  Verificat ion 

It  is  very  difficult  to  measure  the  absorption  of  solar  flux 
within  the  layers  of  the  canopy.  However,  the  model  can  be  bench- 
marked  against  the  measured  reflected  and  transmitted  solar  flux 
densities  within  and  above  the  canopy  system. 

A unique  data  base  exists  for  the  lodgepole  pine  stand  at  Lead- 
ville,  Colorado  (Ranson,  Kirchner,  and  Smith,  1978).  The  data  base 
consists  of  optical  and  thermal  spectral  measurements  for  various 


f 


1 
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terrain  and  temporal  features.  In  this  particular  study,  the  Scene 
Recording  Radiometer  (SRR)  was  principally  utilized  to  obtain  spectral 
reflectance  and  transmittance  measurements  of  the  four  modeling  trees. 

The  SRR  instrument  is  described  by  Berry,  Keir.;es,  and  Smith 
(1978)  . The  SRR  was  suspended  on  support  cables  attached  to  two  15  m 
towers,  which  allowed  spectral  reflectance  measurements  from  above  the 
canopy  to  be  obtained  (Figure  8) . The  SRR  consists  of  a six  narrow 
band  interference  filter  wheel  interfaced  to  a Hasselblad  EL  500 
camera  to  provide  a photographic  record  of  the  scene.  All  filtered 
spectral  data  were  referenced  to  a barium  sulfate  painted  panel  to 
provide  reflectance  values.  Filters  used  were  centered  at  4800,  6750, 

O 

7300,  8000,  and  9600  A.  (The  standard  field  of  view  (FOV)  was  22.5°.) 
Vertical  spectral  reflectance  above  the  four  modeling  trees  was 
measured  throughout  1977  and  compared  with  the  simulated  reflectance 
values . 

On  September  17,  1977  at  1030  Standard  Time  the  SRR  was  placed 
on  the  ground  looking  vertically  up  through  the  canopy  of  the  modeling 
trees.  A 12.7  mm  diameter  stop  on  the  SRR  optics  was  utilized  with 
a FOV  of  9°  which  restricted  the  FOV  to  the  boundaries  of  the  canopy. 
Two  sample  points  were  measured  and  the  mean  transmittance  to  the 
ground  was  calculated. 

Two  experiments  were  performed  to  evaluate  the  spectral  reflec- 
tance variability  of  lodgepole  pine  canopies  with  changing  solar 
zenith  angle.  SRR  data  taken  from  the  tramway  system  of  the  four 
modeling  trees  were  acquired  on  May  24,  1978,  from  0642  - 1526  hours 
Mountain  Daylight  Time  (MDT) . Approximately  10  percent  of  the  scenes 
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a different  cluster  of  three  lodgepole  pine  trees  on  August  4,  1976, 
from  0657  to  1400  hours  MDT . 

Results  and  Discussion 

The  measured  spectral  canopy  transmittance  to  the  ground  versus 
the  simulated  values  are  shown  in  Figure  9.  Considering  the  small 
magnitude  of  the  spectral  transmittance  values,  the.  accuracy  of 
prediction  is  very  good. 

A typical  comparison  between  the  measured  vertical  spectral 

reflectance  of  the  four  lodgepole  pine  trees  versus  the  simulated 

reflectance  is  presented  in  Figure  10. 

The  simulated  total  absorption  coefficients  A of  a single 

z 

isolated  mean  canopy  element  for  solar  zenith  angles  of  0°,  47°,  75°, 

and  85°  were  .68,  .67,  .65,  and  .60,  respectively.  Jarvis  et  al . 

(1976)  states  that  the  optical  properties  of  coniferous  needles  are 

poorly  known  due  to  difficulties  in  making  spectral  measurements  with 

conventional  spectrophotometric  equipment.  However,  Gates  et  al. 

(1965),  assuming  zero  transmissivity,  found  a mean  value  of  absorption 

of  0.88  for  needle  mosaics  of  Pinus  strobus  using  a spectral  solar 

irradiance  curve  of  a sunny  day.  It  is  important  to  note  that  the 

total  A is  not  constant  but  is  a function  of  the  solar  irradiance 
z 

conditions  which  change  as  a function  of  solar  zenith  angle.  Gates 
(1970)  has  shown  that  the  total  absorption  coefficient  of  leaves  of 
various  plants  can  change  as  much  as  0.13  between  sunny  and  cloudy 
irradiance  conditions. 

The  SRVC  absorption  coefficients  cn,  . for  the  six  MCR  for 

A,i,z 

solar  zenith  angles  of  0°,  72°,  and  89°  are  presented  in  Figure  11. 


The  a.  . ° and  a,  . ,°  curves  were  very  similar.  It  is  hypothesized 

X.i.O  X , 1 , 4 7 

that  this  similarity  is  due  to  the  fact  that  the  PGAP  at  solar  zenith 
angles  of  0°  and  45°,  corresponding  to  the  inclination  invervals  of  85 
and  45°,  respectively,  in  Table  4,  are  very  similar  for  any  given 
canopy  layer.  Thus,  the  transfer  of  radiation  should  be  somewhat 
similar.  Several  interesting  trends  are  apparent  in  Figure  11. 

The  trends  seen  for  a.  . can  be  explained  as  follows.  The 
differences  in  the  absorption  coefficients  seen  in  each  canopy  layer 
are  due  to  a complex  interaction  of  radiation,  canopy  geometry, 
optical  properties  of  canopy  elements,  LAI  distributions,  and  spatial 
arrangement  of  canopy  elements.  The  a . ° and  a ^ ° decrease  as 

the  MCR  increases  because  at  high  MCR  an  energy  loss  due  to  high 
element  reflection  is  dominating.  However,  ci,  ° is  relatively 
constant  and  peaks  at  a MCR  of  .3  - .5.  At  the  lower  MCR  relatively 
little  energy  is  reaching  Layer  3 due  to  the  high  element  absorptance, 
and  at  the  higher  MCR  the  high  energy  loss  due  to  high  element  reflec- 
tion is  dominating.  The  ^ q°  (where  4 denotes  the  ground  layer) 
consistently  increases  with  increasing  MCR.  At  higher  MCR  a rela- 
tively large  proportion  of  flux  reaches  the  ground,  and  the  ground 
has  a relatively  high  absorption  coefficient  relative  to  the  canopy 
elements  (Table  3).  These  factors  dominate  as  MCR  increases  and  thus 

a,  , ° increases. 

1,4,0 

It  is  important  that  the  most  drastic  differentials  in  the  a,  . 

X , l , z 

occur  at  low  MCR  which  would  correspond  to  the  photosynthetic  active 
radiation  (PAR)  absorption.  Thus,  the  uneven  distribution  of  absorbed 
PAR  within  the  canopy  will  have  definite  effects  on  the  photosynthetic 
efficiency  of  the  canopy.  In  contrast,  at  high  MCR  the  n.  . 's  tend 

A , 1 , Z 


Table  4.  Simulated  probability  of  gap  (PGAP)  through  each  canopy 

layer  at  the  nine  different  inclination  view  angle  intervals. 


Inclination  view  angle 


5 

15 

25 

35 

45 

55 

65 

75 

85 

Layer  1 

.03 

. 30 

.41 

.46 

.48 

.49 

.50 

.50 

.50 

Layer  2 

.00 

.10 

.19 

.23 

.25 

.26 

.27 

.27 

.27 

Layer  3 

.13 

.47 

.57 

.61 

.62 

.63 

.63 

.63 

.63 
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to  converge  due  to  the  large  degree  of  multiple  scattering,  and 
consequently  a more  even  distribution  of  absorbed  spectral  flux  is 
assumed . 


The  a.  . ,,  curve  increases  slightly  in  magnitude  relative  to 

A , 1 , / Z 

a.  _o  for  all  MCR  which  dictates  that  the  flux  reaching  Layer  2 

A , 1 , U 

should  be  less,  and  as  a consequence  o decreases  slightly  for 

A , Z , / Z 


all  MCR  relative  to  t qo. 


A similar  argument  can  be  made  for  the 


slight  decrease  in  a.  0 -,0o  and  a . 0 which  is  most  noticeable  at 

a j J j / Z A j ‘4  j / Z 

low  MCR.  These  changes  are  not  large,  and  can  be  partially  explained 
by  the  fact  that  the  simulated  PGAP's,  which  are  important  in  radiation 
transfer  for  the  solar  zenith  angles  of  0°,  47°,  and  72°  corresponding 
to  inclination  intervals  of  85,  45,  and  25,  respectively,  do  not 
change  drastically  (Table  4). 

The  ot;  ^ gg0  curve  shows  that  Layer  1 absorbs  most  of  the  energy 
except  at  high  MCR.  At  low  MCR  (e.g.,  PAR)  Layer  1 essentially 
absorbs  all  incident  energy  due  to  the  very  low  PGAP  (Table  4).  The 
a n „ curve  assumes  the  same  general  shape  as  a in  the  previous 

A j ‘•>07  A y J y Z 


Figure  12  presents  the  simulated  proportion  of  spectral  solar 
irradiance  at  various  MCR  values  that  reaches  the  ground  level  at  a 
solar  zenith  angle  of  0°  and  89°.  As  one  proceeds  from  a low  MCR 
to  a high  MCR,  the  proportion  of  spectral  solar  irradiance  reaching 
the  ground  increases  exponentially  due  to  the  increase  in  element 
reflectance  and  multiple  scattering.  This  phenomenon  has  been  observed 
by  many  investigators,  for  example  Jarvis  et  al.  (1976)  and  Ross  (1976). 

Upon  inspection,  the  a in  Figure  11  change  very  little 

A , tOt3l ) Z 

as  a function  of  solar  zenith  angle.  Two  plots  were  produced  for 
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further  analysis.  Figure  13  shows  the  SRVC  simulated  «.  spectral 

1 

absorption  coefficient  for  a MCR  of  0.16  as  a function  of  solar  zenith 
angle.  As  can  be  noted,  the  changes  drastically  as  a function  of 
solar  zenith  angle;  however,  the  total  proportion  of  the  spectral 
irradiance  absorbed  by  the  canopy  system  is  relatively  constant.  The 
modeling  results  from  Chance  and  I.eMaster  (1978)  have  shown  similar 
trends  for  a wheat  canopy.  Figure  13  shows  the  SRVC  simulated 
spectral  absorption  coefficient  for  a MCR  of  0.64  as  a function  of 
solar  zenith  angle.  a.  does  not  change  drastically  with  solar  zenith 
angle  due  to  the  high  element  reflectance  and  a relatively  high  degree 
of  multiple  scattering.  Again  the  total  proportion  of  the  spectral 
irradiance  absorbed  by  the  canopy  system  is  relatively  constant. 

Since  the  total  proportion  of  spectral  irradiance  absorbed  by 
the  canopy  system  changes  little  as  a function  of  solar  zenith  angle, 
the  total  reflected  spectral  flux  would  also  change  little  since  the 
canopy  system's  reflectance  and  absorptance  must  equal  1.0.  However, 
the  SRVC  model  simulates  the  apparent  reflectance  in  the  nine  incli- 
nation bands  above  the  canopy,  and  this  reflectance  is  by  no  means 
constant  with  varying  solar  zenith  angle  due  to  the  complex  radiation- 
canopy  geometry  interactions.  Many  studies  do  not  distinguish  between 


total  and  directional  reflectance  factors  which  may  dictate  what  kind 
of  reflectance  versus  solar  zenith  angle  trends  result  in  any  partic- 
ular study. 

For  example,  within  this  study  the  simulated  total  spectral  flux 
absorbed  by  the  canopy  system  is  essentially  constant  with  solar 
zenith  angle.  However,  the  simulated  vertical  spectral  reflectance 
shows  definite  trends  as  seen  in  Figure  14.  Sun  angle  experiments. 


4 


Figure  13.  SRVC  simulated  spectral  absorption  coefficients  for  a 
mean  canopy  component  reflectance  (MCR)  of  0.64  and  0.16 


of  three  trees,  and  the  measured  data  of  the 
modeling  trees. 
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discussed  pi'eviously,  using  the  SRR  instrument  were  conducted  on  the 
four  modeling  trees  and  on  a nearby  cluster  of  three  trees.  The 
results  of  vertical  reflectance  versus  solar  zenith  angle  of  these 
two  targets  for  the  0.68  pm  band  and  the  0.80  pm  band,  corresponding 
to  a MCR  of  .064  and  .53S,  respectively,  are  presented  in  Figure  14. 

The  simulated  results  of  the  three  trees  deviate  from  the  measured 
results  in  magnitude,  which  might  be  expected  since  the  targets  are 
two  distinct  points.  However,  a decreasing  reflectance  with  increasing 
solar  zenith  angle  is  seen  in  both  simulated  and  measured  results. 

The  same  general  trends  exist  for  the  results  of  the  four  modeling 
trees  with  better  correspondence  in  magnitude  to  the  simulated 
results  (Figure  14) . 

The  simulated  proportion  of  total  global  solar  irradiance  absorbed 
by  the  lodgepole  pine  canopy  system  as  a function  of  solar  zenith 
angle  is  presented  in  Figure  15.  The  total  global  solar  irradiance 
absorbed  by  the  canopy  system  is  essentially  constant  with  changing 
solar  zenith  angle,  and  therefore,  the  total  reflectance  is  essentially 
constant.  However,  the  proportions  absorbed  by  Layers  1 and  2 are 
clearly  dependent  on  solar  zenith  angle.  As  solar  zenith  angle 
increases,  the  PGAP  for  Layer  1 decreases  and  thus,  the  absorption  in 
Layer  1 increases.  In  contrast,  as  solar  zenith  angle  increases,  the 
proportion  of  flux  reaching  Layer  2 decreases  due  to  the  lower  PGAP 
of  Layer  1 at  higher  solar  zenith  angles. 

Bergen  (1971,  1974)  has  shown  in  a 10  m tall  lodgepole  pine 
canopy  that  the  air  temperature  maximum  occurs  in  the  upper  crown 
during  the  morning  period  and  then  descends  deeper  into  the  crown  as 
the  solar  zenith  angle  decreases  for  clear  sky  and  low  wind  speed 


Figure  15.  Simulated  proportion  of  global  solar  irradiance  absorbed 
by  the  lodgepole  pine  canopy  system  (total).  Layer  1, 
Layer  2,  Layer  3,  and  the  ground,  as  a function  of  solar 
zenith  angle  for  October  14,  1977. 
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conditions.  The  maximum  temperatures  move  from  the  9 m level  to  the 
7 m level  during  the  day.  This  phenomenon  suggests  that  the  foliage 
heated  by  the  solar  irradiance  is  the  cause  of  the  shifting  maximum 
air  temperature  level  (Bergen,  1971).  Figure  15  demonstrated  a shift 
of  levels  at  which  the  maximum  global  solar  irradiance  is  absorbed. 

At  the  greatest  solar  zenith  angles,  Layer  1 clearly  absorbs  the 
largest  proportion  of  global  solar  irradiance.  However,  at  relatively 
low  solar  zenith  angles,  the  proportion  of  solar  global  irradiance 
absorbed  by  Layers  1 and  2 begins  to  converge  indicating  the  level  of 
maximum  absorbed  solar  irradiance  is  shifted  down  into  the  canopy. 

This  can  be  explained  by  the  PGAP  as  a function  of  view  angle  as 
discussed  above.  Bergen  (1974)  states  that  forest  productivity 
studies  require  more  quantitative  data  of  canopies  such  as  radiation 
flux  absorption. 

In  future  work  it  is  anticipated  that  the  SRVC  absorption  model 
will  be  used  to  study  PAR.  absorption  within  canopies  of  various 
geometric  structures.  In  addition,  the  results  of  this  study  will  be 
combined  with  a comprehensive  thermal  canopy  signature  model  (TCSM) 
recently  developed  at  Colorado  State  University  to  study  the  thermal 
behavior  of  the  lodgepole  pine  canopy. 

Conclusions 

A solar  radiation  vegetation  canopy  absorption  model  has  been 
developed  which  has  several  advantageous  features  for  quantifying 
solar  absorption  within  vegetation  canopies  under  a variety  of  environ- 
mental and  plant  conditions.  Under  a theoretical  clear  and  dry 
atmosphere  and  a correspondingly  low  diffuse/direct  solar  irradiance 


ratio,  the  following  measured  and  simulated  results  were  found  for  a 

specific  lodgepole  pine  canopy. 

1)  The  simulated  results  correspond  relatively  well  with  the 
measured  reflectance  and  transmittance  data. 

2)  Large  differentials  occur  in  spectral  absorption  by  canopy 
layers  especially  in  the  photosynthetic  active  radiation  region 
as  a function  of  solar  zenith  angle  and  canopy  geometry,  which 
may  have  significant  effects  on  the  photosynthetic  efficiency 
of  the  canopy. 

3)  As  one  proceeds  from  a low  MCR  to  a high  MCR  the  proportion 
of  spectral  irradiance  reaching  the  ground  increases 
exponentially. 

4)  Although  it  is  believed  that  the  total  spectral  reflection  by 
the  canopy  system  is  relatively  constant  with  solar  zenith 
angle,  the  vertical  spectral  reflectance  of  the  canopy 
decreases  with  increasing  solar  zenith  angle. 

5)  Both  the  total  and  spectral  global  irradiance  absorption  by 
the  entire  lodgepole  pine  canopy  system  are  relatively  constant 
with  solar  zenith  angle.  However,  the  proportion  of  total  and 
spectral  global  irradiance  absorbed  by  individual  canopy  layers 
varies  greatly  as  a function  of  solar  zenith  angle. 
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3.0  A THERMAL  VEGETATION  CANOPY  MODEL  OF  SENSOR  RESPONSE 


Abstract 

A thermal  canopy  signature  model  (TCSM)  was  developed  to  approx- 
imate the  thermal  behavior  of  a vegetation  canopy  by  a mathematical 
abstraction  of  three  horizontal  layers  of  vegetation.  The  geometry 
of  canopy  elements  within  each  layer  is  quantitatively  described  by 
the  foliage  orientation  distribution  and  number  density.  Given  this 
geometric  information  for  each  layer  and  the  driving  variables  (direct/ 
diffuse  solar  irradiatice,  air  temperature,  horizontal  wind  velocity, 
relative  humidity,  and  ground  temperature)  the  energy  budgets  of  average 
leaves  within  each  layer  are  determined.  The  resulting  system  of 
conservation  equations  is  solved  for  the  average  layer  temperature. 

This  information  is  then  used  to  calculate  the  response  of  a thermal 
infrared  sensor  at  varying  view  angles  above  and  within  the  canopy. 

The  model  is  applied  to  lodgepole  pine  (Pinus  contorta)  canopy  and 
the  results  are  validated  with  both  radiometric  and  contact  temper- 
ature measurements.  The  simulated  average  layer  temperatures 
closely  follow  air  temperature  due  to  the  high  leaf  area  index  values 
of  the  lodgepole  pine  canopy  and  the  small  dimensions  of  the  needles. 
Simulated  versus  measured  radiometric  average  temperatures  of  Layer  2 
correspond  approximately  within  2°C.  Simulated  results  suggest  that 
canopy  element  geometry  can  significantly  influence  the  effective 
radiant  temperature  of  a sensor  above  the  canopy  as  a function  of  view 
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angle.  This  phenomenon  has  important  implications  on  the  optimum 
view  angle  for  making  inferences  about  the  target  of  interest. 

Introduction 

Radiant  energy  interacts  with  a vegetation  canopy  in  a complex 
manner.  A vegetation  canopy  can  generally  be  defined  as  a stand  of  one 
or  many  plant  species.  The  interactions  of  radiation  and  other  energy 
exchanges  with  the  individual  canopy  elements  (leaves,  stems,  and 
reproductive  structures)  include  solar  radiant  absorption,  reflection, 
and  transmission;  thermal  radiant  emission,  absorption,  and  reflec- 
tion; convection;  transpiration;  conduction;  photosynthesis;  and 
respiration.  The  interactions  between  the  ensemble  of  canopy  elements, 
the  sky,  ground,  and  air  determines  the  resulting  thermal  behavior, 
solar  radiation  regime,  and  thermal  radiation  regime  of  the  canopy 
structure.  Ross  (1976)  discusses  the  factors  which  determine  the 
resulting  radiation  regimes  in  plant  canopies.  Monteith  (1973)  and 
Gates  (1963)  discuss  the  thermal  behavior  of  vegetation  elements. 

An  understanding  of  the  underlying  principles  involved  in  the 
above  radiation  and  energy  regimes  is  important  to  remote  sensing 
applications.  Such  basic  knowledge  is  needed  to  improve  the  accuracy 
of  remote  sensing  techniques  for  determining  species  identification, 
vegetation  stress,  vegetation  biomass  and  the  nature  of  the  underlying 
soil  and/or  rock. 

One  method  used  to  study  these  interactions  is  by  mathematical 
modeling  techniques.  Canopy  modeling  enables  the  experimentalist  to 
conveniently  organize,  in  a mathematical  sense,  all  the  complex  inter- 
actions which  take  place  in  a vegetation  canopy,  and  thus  enables  one 
to  integrate  information  and  interpret  the  data  collected.  In 
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addition,  models  can  serve  as  a guide  to  experimentation.  As  in  any 
scientific  research,  knowledge  of  the  underlying  basic  principles  of 
a phenomenon  needs  to  be  understood  before  a surge  of  practical 
applications  can  be  realized. 

Several  solar  radiation  canopy  models  which  incorporate  canopy 
geometry  exist.  Some  of  the  more  recent  developments  have  been  made 
by  Oliver  and  Smith  (1974),  Suits  (1972),  Idso  and  deWit  (1970),  and 
Allen  and  Richardson  (196S).  It  is  known  that  the  spectral  signature 
of  most  vegetation  canopies  varies  with  both  direction  of  sensor  view 
and  solar  zenith  angle.  This  variation  is  primarily  due  to  differences 
in  canopy  geometry  which  influence  the  transfer  of  radiation  within 
a vegetation  canopy.  The  canopy  geometry  can  be  described  by  such 
physical  characteristics  as  the  distribution  of  plants  on  the  ground, 
leaf  area  index  and  its  distribution  as  a function  of  height,  leaf 
angle  frequency  distribution  and  leaf  azimuth  angle  frequency  distri- 
butions. These  canopy  characteristics,  in  regard  to  radiation 
transfer,  are  discussed  by  Oliver  and  Smith  (1974),  Idso  and  deWit 
(1970),  and  deWit  (1965). 

The  thermal  infrared  region  (3-20  pm)  of  the  electromagnetic 
spectrum  in  recent  years  has  received  keen  interest.  This  radiation 
region  may  add  valuable  additional  information  to  make  inferences 
concerning  the  characteristics  of  vegetation  canopies.  Many  thermal 
models  exist  for  different  non-vegetated  targets  of  interest.  Several 
models  exist  for  planar  solid  objects.  For  example,  Watson  (1971) 
developed  a thermal  model  for  predicting  the  diurnal  surface  tempera- 
ture variation  of  the  ground,  and  the  University  of  Michigan  (1969) 
developed  a model  for  the  prediction  of  time-dependent  temperatures 


and  radiance  of  planar  targets  and  backgrounds.  However,  few  thermal 
models  exist  for  plant  canopies. 

Gates  (1968)  presents  an  energy  budget  for  a single  plant  leaf 
isolated  in  space.  In  addition,  Rimes,  Ranson,  Kirchner,  and  Smith 
(1978),  and  Wiebelt  and  Henderson  (1977)  have  developed  thermal  models 
of  an  individual  leaf.  Other  investigators  model  the  thermal  dynamics 
of  vegetation  canopies  assuming  a simplistic  single  homogeneous  layer 
abstraction.  For  example,  vegetation  is  treated  as  a single  homo- 
geneous layer  with  an  associated  transmission  factor  for  solar  radia- 
tion in  the  University  of  Michigan  model  (1969)  . Heilman  e_t  a_l.  (1976) 
used  actual  thermal  scanning  data  to  measure  crop  effective  radiant 
temperatures  and  used  an  evapotranspiration  (ET)  equation  to  estimate 
crop  ET.  However,  they  assume  they  are  viewing  only  the  top  layer  of 
the  crop  with  the  scanner  and  ignore  the  canopy  geometry.  The  litera- 
ture review  failed  to  reveal  any  thermal  canopy  models  which  physically 
account  for  the  canopy  geometry  and  the  thermal  dynamics  within  the 
canopy. 

It  is  known  that  vegetation  canopies  are  non-Lamber tian  at  optical 
wavelengths,  primarily  due  to  canopy  geometry.  Similarly,  in  the 
thermal  region,  it  is  believed  that  while  individual  canopy  elements 
are  isotropic  radiators,  the  response  from  the  canopy  may  also  be 
non-bambertian  because  canopy  geometry  causes  spatial  variations  in 
many  energy  flow  processes. 

With  the  number  of  thermal  sensor  systems  currently  operating  in 
present  satellites  (e.g..  Heat  Capacity  Mapping  Mission)  and  proposed 
for  future  satellites  (e.g..  Thematic  Mapper  on  Landsat  D) , thermal 
models  will  become  increasingly  important  in  interpreting  the  resulting 
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data.  Vegetation  is  often  the  target  of  interest,  especially  in 
agriculture  and  forestry  applications.  However,  in  many  cases  the 
substrate  (e.g.,  soil,  rock)  underlying  the  vegetation  is  of  interest. 

Thermal  modeling  of  vegetation  canopies  could  play  an  important  role 
in  interpreting  thermal  data  (Watson,  Rowan,  and  Of  field,  1971)  and 
in  design  studies  for  discriminating  between  vegetation  types  or 
background  materials  under  a variety  of  environmental  conditions. 

Mo del  Description 

The  primary  objectives  of  this  study  were  to  produce  a thermal 
canopy  model  which  simulates,  in  a physically  based  manner,  (1)  the 
geometric  arrangement  of  primary  canopy  elements,  (2)  the  decreased 
direct/dif f use  solar  radiation  absorption  due  to  the  scattering  of 
neighboring  canopy  components,  (3)  the  increased  thermal  absorption 
of  leaves  due  to  the  thermal  emissions  of  neighboring  canopy  components, 

(4)  the  true  average  temperature  of  scene  elements  within  three  hori- 
zontal, infinite  canopy  layers,  and  (5)  the  response  of  a thermal 
sensor  at  varying  view  angles  above  the  canopy  and  at  horizontal 
looking  positions  within  the  canopy. 

Canopy  Abstraction 

The  vegetation  canopy  is  abstracted  as  three,  statistically 
independent,  horizontal,  infinite  layers  (Figure  1).  The  canopy 
elements  (e.g.,  leaves,  branches,  and  other  plant  organs)  within  each 
layer  are  described  as  a statistical  ensemble  which  is  used  to  define 
the  canopy  geometry.  Mid-elements  that  represent  canopy  elements 
which  occur  at  the  horizontal  plane  occupying  the  middle  of  each 
layer  are  defined.  An  energy  budget  equation  is  formulated  for  the 
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mid-elements  of  each  layer.  These  equations  account  for  the  energy 

inflow  and  outflow  processes  of  the  elements.  The  energy  transfers 

_2 

are  calculated  on  a power  per  unit  area  of  element  (wm  “)  basis. 

The  roots  of  the  resulting  system  of  equations  are  the  average  surface 
temperature  of  the  mid-elements  in  the  three  layers.  It  is  assumed 
that  these  values  represent  the  average  temperature  of  the  elements  in 
each  respective  layer.  These  values  can  then  be  utilized  to  calculate 
the  response  of  a thermal  sensor  at  varying  view  angles. 

The  flow  of  energy  within  a canopy  is  time  dependent.  However, 
the  TCSM  assumes  a steady-state  condition  in  which  elements  of  the 
canopy  are  neither  gaining  nor  losing  a net  amount  of  energy.  In 
addition,  the  energy  loss  due  to  photosynthesis  and  energy  gain  by 
respiration  is  assumed  negligible  and  has  been  ignored.  Heat  exchange 
by  conduction  is  also  considered  negligible,  and  all  surfaces  of 
finite  elements  within  a layer  are  considered  to  be  the  same  tempera- 
ture. These  approximations  are  good  for  elements  of  relatively  small 
dimensions  (Gates,  1975).  The  above  steady-state  and  conduct  on 
assumptions  may  not  he  adequate  when  dealing  with  canopies  which 
exhibit  a large  fraction  of  the  total  element  surface  areas  as  rela- 
tively large  branches  and  trunks.  To  approximate  time  dependent 
events,  one  can  consider  a series  of  incremental  changes  in  steady- 
state  energy  flow  as  discussed  later. 

Several  other  assumptions  are  made.  First,  the  spectral  effects 
in  the  thermal  region  are  assumed  insignificant.  Kondrat'ycv  (1965) 
states  that  natural  surfaces  can  he  treated  in  the  first  approximation 
as  gray  body  radiators  and  emitters.  Data  from  Lecman  et  al^.  (1971) 

I 
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show  that  the  thermal  spectral  emlssivity  of  plants  are  essentially 
constant  with  wavelength. 

Secondly,  the  reflection  of  thermal  flux  within  the  canopy  is 
ignored.  Ross  ( 1976)  states  that  the  transfer  theory  for  the  thermal 
radiation  in  a vegetation  canopy  differs  from  shortwave  theory  in  that 
the  scattering  of  thermal  radiation  may  be  neglected  but  the  emission 
of  thermal  radiation  from  plant  elements  must  be  acknowledged.  It  is 
believed  that  within  natural  vegetation  canopies,  reflected  thermal 
radiation  is  a negligible  contribution  to  the  total  energy  budget. 
Blaxter  (1967)  reported  the  enissivity  of  "green  grass"  as  0.99. 

Idso,  Baker  and  Blad  (1969)  reported  the  emissivity  of  34  plants 
ranging  approximately  from  0.94  to  1.00,  with  30  of  the  plants  above 
0.96.  The  effect  cf  ignoring  thermal  reflectance  on  the  final  temp- 
erature of  a single  leaf  was  explored  using  the  following  analogy 
(Figure  2).  The  walls  of  the  box  are  considered  to  be  leaves  which 
have  an  emissivity  of  .95.  The  energy  budget  is  used  to  calculate  the 
leaf  temperature  of  the  theoretical  box  in  which  thermal  reflections 
are  complete  both  for  the  exterior  and  interior  of  the  box,  simulating 
a theoretical  enclosed  canopy.  Two  calculations,  one  ignoring  and  one 
accounting  for  thermal  reflectance,  are  made  of  the  single  leaf  temp- 
erature which  is  completely  enclosed  in  the  box  composed  of  other 
leaves.  It  is  believed  the  difference  between  these  two  calculations 
will  simulate  the  worst  possible  case  for  ignoring  thermal  reflectance. 
Within  a normal  canopy  a single  leaf  is  rarely  completely  surrounded 
by  neighboring  leaves  and  thus  some  of  the  reflected  thermal  radiation 
escapes  out  the  sides  of  the  system.  However,  this  is  not  the  case  in 
the  "black  needle  box."  The  energy  budget  with  the  specific 
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coefficients  used  to  calculate  the  above  leaf  temperature  is  presented 

in  detail  by  Rimes,  Ranson,  Kirchner,  and  Smith  (1978). 

Several  permutations  of  environmental  conditions  were  simulated 

for  this  simple  example.  The  constants  were:  wind  velocity  1 cm/sec, 

leaf  dimension  1.0  cm,  internal  leaf  resistance  to  water  vapor 

diffusion  .66  min/cm,  leaf  emissivity  0.95.  Simulations  were  run  with 

permutations  of  air  temperature  (0,  15,  35°C),  relative  humidity  (10, 

50,  100%) , accounting  for  and  ignoring  transpiration,  and  no  solar 

2 

irradiance  and  336  w/m  absorbed  solar  irradiance  by  the  leaf.  The 
results  showed  that  the  maximum  difference  between  the  two  calculations 
for  all  permutations  was  1.59°C  ±0.01. 

Finally,  the  individual  canopy  elements  are  assumed  to  emit 
thermal  radiation  in  an  isotropic  manner.  Kondrat'yev  (1965)  and 
Hudson  (1969)  state  that  the  radiation  emitted  from  natural  surfaces 
is  essentially  isotropic. 

In  the  following  discussion  a description  of  the  canopy  geometry 
will  be  presented  followed  by  the  energy  budget  equations  for  each 
layer  which  account  for  the  thermal  radiation  transfers,  solar  radia- 
tion absorption,  thermal  cxitance,  transpiration  and  convection 
exchanges . 

Canopy  Geometry 

Important  parameters  in  describing  radiation  transfer  in  complex 
structures  are  the  gap  frequency  and  the  extinction  of  radiation  with- 
in the  structure.  Monteith  (1965),  Warren  Wilson  (1965),  deWit  (1965), 
and  other  authors  have  developed  various  formulas  for  these  parameters. 


I 


3-10 


Nilson  (1971)  presented  a good  review  of  these  formulations  for 
theoretical  models  of  canopy  geometry  which  have  been  utilized. 

The  geometry  of  the  thermal  canopy  model  is  abstracted  in  the 
following  manner.  Since  the  model  is  numerical  as  opposed  to  analyt- 
ical in  nature,  the  hemispheres  above  and  below  a particular  layer 
are  discretized  into  9 hemispherical  inclination  bands  0-90  degrees 
(Figure  3).  Each  of  the  9 bends  is  further  discretized  into  18 
azimuthal  sectors  (Figure  3).  Within  each  sector  the  radiation  trans- 
fers between  the  three  canopy  layers,  ground  and  sky  are  calculated. 

The  formulation  developed  by  Idso  and  deWit  (V  0)  has  been 
incorporated  to  predict  the  probability  of  gap  in  the  direction  of 
the  nine  hemispherical  bands  for  each  of  the  three  layers. 

The  positive  binomial  distribution  is  used  to  descr  1 .=>'  probabil- 

ities and  azimuthal  symmetry  is  assumed.  The  probab.iity  of  gap  in  a 
particular  band  direction  is  equal  to  the  ratio  of  the  projection  of 
elements  in  a layer  to  the  projection  of  the  underlying  soil  surfaces. 
For  a hemispherical  band  direction  j the  equation  is: 


PGAP . 


jm 
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S -G . 
m i rn 
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m ; PH1T.  = 1 - PGAP. 

jm  jm 


where : 


PGAP j ^ = probability  of  gap  for  layer  m in  direction  of 
hemispherical  band  j 

PHIT.  = probability  of  hit  for  layer  in  in  direction  of 
jm 

hemispherical  band  j 

G.  = mean  canopy  projection  of  elements  in  layer  m in 

the  direction  of  hemispherical  band  j 
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Horizontal  view 


SOUTH 


V l ow 


Horizontal  and  vertical  views  of  the  9 hemispherical 
inclination  hands  which  are  divided  into  18  equal 
azimuthal  sectors.  Tf  one  rotates  the  horizontal  view 
about  the  7.  axis  the  hands  would  occur  in  three-snace . 
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l.A  1 = element  area  index  for  layer  m 

m 

= index  of  spatial  dispersion  of  elements  in  layer  m 
sin(0.)  = sine  function  of  the  inclination  angle  0 of 
hemispherical  band  j. 

The  function  G.  is  determined  from  inclination  angle  frequency  dis- 
j r.i 

trihutions  of  the  elements  in  a layer.  The  derivation  and  computa- 
tional procedure  is  presented  by  deWit  (1965).  The  parameter 
ranges  from  0 to  1 and  is  an  index  of  denseness  or  spatial  dispersion 
of  the  el  ments  in  a canopy.  As  S approaches  1.0,  the  more  regular 
the  dispersion  of  elements  is  and  the  less  frequently  a gap  is 
encountered.  The  leaf  area  index  of  a canopy  layer  is  equal  to  the 
ratio  of  the  total  one-sided  element  area  within  a layer  to  the  area 
of  the  underlying  soil  area.  For  a more  in-depth  discussion  of  the 
above  theory  and  the  required  measurements  see  Idso  and  deWit  (1970), 
and  Oliver  and  Smith  (1974). 

The  resulting  PGAP  ^ ^ and  PllIT.^  are  important  parameters  in 
describing  the  radiation  transfers  with  each  hemispherical  sector. 

In  addition,  the  probabilities  of  gap  and  hit  of  half  of  each  layer 
are  required,  and  these  parameters  are  calculated  as: 


PGAP = (PGAP  . ) 
dm  jm 


PH  IT', 
jm 


1 - PGAP! 

jm 


where: 


PGAPj^  = probability  of  gap  for  one  half  of  layer  m in 
the  direction  of  hemispherical  band  j 
PH IT l ^ = probability  of  hit  for  one  half  of  layer  m in 
the  direction  of  hemispherical  band  j. 


t 
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Thermal  Radiation  Transfers 

The  following  describes  the  manner  in  which  thermal  radiation 

transfers  are  calculated  within  the  model.  Each  layer  emits  and 

receives  thermal  radiation  in  the  hemispheres  occurring  above  and 

below  a particular  layer.  The  transfer  of  thermal  radiation  within 

each  hemispherical  sector  between  the  three  canopy  layers,  the  sky, 

and  ground  is  calculated  as  follows.  As  seen  in  Figure  4,  for  small 

angles  (0,<J0  the  two  sides  of  a sector  can  be  described  as  r cos0  d 4> 

and  r d0,  and  the  area  of  the  sector  is  described  as  approximately 
2 

r cosO  d0  d<}>.  One  can  then  define  the  solid  angle  of  a sector  as: 

ho  - r^  c.osO  dO  df 

2 

r 

= cosO  d0  d<{> 

where : 


£1  = steradians  of  a sector. 


And  it  follows  that 

-0. 


ft 


whe  re : 


-/  / 
J0. 


cosO  dG  d'j> 


<t>l»  <J>2  = define  the  azimuthal  limits  of  sector  i in 
hemispherical  band  j 

0^,  02  = define  the  inclination  limits  of  sector  i in 
hemispherical  band  j . 

To  calculate  the  thermal  irradiance  on  a mid-leaf  from  a parti- 
cular layer  in  any  given  sector  we  proceed  as  follows.  Assuming  that 
canopy  elements  in  a particular  layer  emit  thermal  radiation  in  an 
isotropic  manner,  the  radiance  from  the  material  is 
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4 


L = — 

7T 


(i) 


where: 


-2  -1 

L = radiance  (wm  *sr  ) 


M = exitance  (wm  “). 

The  above  radiance  L is  equal  for  all  viewing  directions;  however,  a 
canopy  layer  has  special  characteristics  in  that  it  is  not  solid  but 
has  gaps  which  are  dependent  on  the  direction  of  view.  As  a conse- 
quence, the  irradiance  on  a panel  normal  to  the  mid-vector  (Figure  A) 
of  a sector,  defined  by  0^,  and  and  ^ from  an  infinite  hori- 
zontal layer,  is  calculated  by 


rr 

1J  Ja.  h. 


where: 


L*  PHIT (0)  cosO  dQ  dij> 


-2s 


L 


= the  irradiance  (wm  ) on  a panel  normal  to  the 
mid-vector  from  the  sector  i in  hemispherical 
band  j 

= radiance  of  canopy  elements 
PHIT(O)  = probability  of  hit  for  viewing  angle  0. 

The  above  assumes  that  elements  within  a layer  have  a homogeneous 
surface  temperature  and  emissivity.  The  equations  and  theory  of  flux 
transfer  from  extended  sources  through  solid  angles  to  receiving 
elements  are  presented  by  the  National  Bureau  of  Standards  (1978). 


Assuming  that  PHIT(0)  is  constant  within  sector  ij , then 

cosO  dO  d<j>. 


*2  r°2 


E.  . = L-PHIT 
1J 


.rr 

JA.  Je. 


A *■  v •»» 
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Because  the  18  sectors  within  band  j have  equal  solid  angles,  the 
above  equation  can  be  reduced  to: 

L PH  IT , r 2;t  r0n 


E.  . = 
1 1 


18 


77' 

Jo  Jo, 


cos9  d9  d$ 


■PH IT. 
1 


ij 


if 

J 0, 


cosO  d0. 


This  expression  can  be  further  evaluated  as: 

(sin9?  - sin9  ) 

E..  = L’PHIT . * it  — (2) 

ij  j 9 

Combining  equations  (1)  and  (2)  for  a particular  sector 

(sin9„  - sinO  ) 

E = M-PHIT  . — - 

3 o 


For  simplicity,  the  quantity  (sinC^  - sinO  )/9  will  be  defined  by 
SECTOR j where  j denotes  the  hemispherical  band  interval: 

E..  = M*PHIT . • SECTOR, . 
ij  J J 

The  above  assumes  that  the  panel  which  represents  a mid-elen.ent 
is  normal  to  the  direction  of  the  source  and  that  there  exist  no 
obstructions  between  the  emitting  canopy  layer  and  the  panel.  The 
following  calculations  correct  for  the  fact  that  the  panel  or  mid- 
element is  not  always  oriented  normal  to  the  source. 

The  desired  correction  factor  is  the  cosine  of  the  angle  between 
the  source  vector  and  the  normal  vector  of  the  panel.  The  theory  is 
based  on  the  existence  of  planar  elements.  The  inclination  angles  of 
the  canopy  elements  and  source,  and  the  azimuthal  angles  of  the  leaves 
and  source  are  discretized  as  before.  The  canopy  elements  are  assumed 
to  have  azimuthal  symmetry.  The  direction  cosines  of  all  source 


sectors  are  calculated  as: 


xs. . 

1J 

YS.  . 

= 

XS.  . 

ij 

si"<V 


-.share: 

Vg  = vector  of  direction  cosines  for  source  sector  i 

in  hemispherical  band  j 

0. . and  <p. . = the  inclination  and  azimuth  angle,  respectively, 
ij  ij  7 ’ 

of  the  mid-vector  in  sector  i and  hemispherical 
band  j . 

The  direction  cosines  for  the  normal  vector  of  all  planar  element 
inclination  angle  intervals  are  calculated  as  follows.  The  azimuth 
angle  is  fixed  to  zero  degrees  since  the  canopy  is  assumed  azimuthally 
symmetric,  both  in  geometric  and  thermal  radiant  energy  modes.  Thus, 
regardless  of  the  azimuthal  orientation  of  an  clement,  the  thermal 
radiant  contributions  to  the  element  are  constant  for  any  specific 
element  inclination  angle. 

XL^  -sin(O^) 

V = YL.  = 0.0 

L J 

ZLj  cos(Oj) 

where: 

= vector  of  direction  cosines  for  the  normal  of  a planar 
element  with  inclination  j 

0^  = element  inclination  of  hemispherical  band  j. 

Now  one  can  calculate  the  absolute  value  of  the  dot  products  for 
all  source-element  angle  permutations.  These  values  are  equal  to  the 


a k.  ’ll 
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correction  factors  desired. 

c"5uk  ■ 'V\l 

where: 


COS  = the  correction  factors  desired  for  permutations 

ijk 

of  source  sector  i in  hemispherical  band  j and 
element  inclination  k. 

In  addition,  one  must  apply  the  absorption  coefficient  for  thermal 
radiation.  As  a result,  the  equation  becomes 


<t>.  4 1, 

i .lk  _ 


m 


M-PHIT. • SECTOR. -ABSORB -COS. 

J J i J k 


where: 

$ , 

— = tlie  thermal  flux  density  (vm  ) absorbed  by  a 
m 

mid-element  inclined  at  inclination  angle  k from 
source  sector  i in  hemispherical  band  j 
ABSORB  = thermal  absorption  coefficient. 

The  above  assumes  that  there  exist  a single  layer  and  a removed 
single  element  receiving  flux  from  that  layer.  However,  the  contri- 
bution of  absorbed  thermal  flux  density  from  all  hemispherical  sectors 
both  upward  and  downward  directions  for  each  canopy  layer,  the  sky, 
and  ground,  to  each  layer's  mid-elements  must  be  calculated  (Figure  5) 
The  calculations  should  account  for  the  fact  that  within  each 
sector  the  flux  which  originates  from  any  given  layer  is  obstructed 
by  other  leaves  before  the  flux  reaches  any  specified  mid-element  in 
another  layer.  In  addition,  a relatively  large  number  of  permutations 
must  be  calculated  since  each  layer  is  simultaneously  emitting  thermal 
flux  to  other  layers  and  absorbing  emitted  flux  from  the  surrounding 
leaves,  other  layers,  the  sky,  and  the  ground.  For  each  permutation 


Figure  5.  Hemispherical  sectors  are  shown  for  the  sky. 
Layer  1,  and  the  ground.  Note  only  sectors 
in  one  band  are  shown. 
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a contribution  coefficient  which  replaces  HUT.  is  calculated.  For 

J 

example,  the  mid-elements  in  Layer  1 will  receive  thermal  flux  from 
the  sky.  Layer  1,  Layer  2,  Layer  3,  and  the  ground.  For  all  sectors 
defined  within  a specific  hemispherical  band  j,  the  contributing 
coefficients  are  calculated  as  follows. 

The  proportion  of  sky  thermal  flux  within  a sector  in  band  j 
reaching  the  mid-elements  in  Layer  1 is 


PGAP'  . 

Jl 

The  contributing  coefficient  from  Layer  1 to  the  Layer  1 


mid-elements  is 


2*(PKIT') 


The  coefficient  of  2 represents  the  two  half-layers  of  Layer  1. 

The  contributing  coefficient  from  Layer  2 to  the  Layer  1 mid- 
elements can  be  derived  in  the  following  manner.  The  probability  of 
gap  to  Layer  2 is  PGA?',  j.  Once  Layer  2 is  reached,  the  projected 
surface  area  of  interest  is  PlirT.?.  Thus,  the  contributing  coeffi- 


cient is 


PGAP , *P1I IT... 
Jl  J2 


A similar  argument  can  be  made  for  the  contributing  coefficient 

from  Layer  3 to  the  Layer  1 mid-elements: 

PGAP'  • PGAP.. -PH IT... 

J 1 J 2 .1 3 

The  contributing  coefficient  from  the  ground  to  Layer  1 mid- 
elements is 

PGAP  ’. . • PGAP  . „ • PGAP  . , . 

J 1 .1  2 J 3 
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The  contributing  coefficients  for  both  upward  and  downward  direc- 
tions of  a particular  sector  and  layer  should  sum  to  2.0  representing 
the  two  sides  of  the  mid-elements.  And,  in  fact,  if  one  sums  the 
above  coefficients,  the  total  is  2.0. 

In  a similar  fashion,  the  contributing  coefficients  from  all 
source  sectors  to  Layers  2 and  3 mid-elements  are  calculated. 

The  final  equation  which  calculates,  within  a particular  sector, 
the  amount  of  flux  density  absorbed  by  a mid-element  at  a particular 
inclination  angle  from  any  given  source  layer  is 
$ . 

.JJiiEl  = >(  -coNT . ’SECTOR.  ‘ABSORB  -COS.., 

2 i j im  ] m ljk 

m j j j 

where : 


v. 

ij  k m 

2 

m" 


= within  source  sector  i in  hemispherical  band  j,  it 
is  the  thermal  flux  density  absorbed  by  a mid- 
element in  layer  m inclined  at  inclination  angle  k 
from  source  elements  in  layer  i.  Note  the  index  i 
represents  the  sky  and  ground  in  addition  to  the  3 


canopy  layers. 


M = average  thermal  exitance  of  elements  in  layer  i 

CONTj ^ = contributing  coefficient  for  mid-elements  in  layer  m 
absorbing  flux  from  elements  in  layer  i for  all 
sectors  within  hemispher ical  band  j 
ABSORB^  = average  thermal  absorption  coefficient  for  elements 
i n 1 aye  r tn . 

Except  for  the  sky  thermal  exitance,  M can  be  further  expressed  in 
terms  of  the  Stef an- Boltzmann  Law: 

M = o«e  *T  4 

l l l 
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where: 

0 = Stefan-Boltzmann  constant 

= average  emissivity  of  elements  in  layer  \ 

T = true  average  surface  temperature  (°K)  of  the  mid-elements 
in  layer  x (unknown) . 

Note  that  the  average  surface  temperatures  of  the  three  canopy  mid- 
elements are  not  known  and  these  values  must  be  derived  mathematically. 
In  addition,  the  ground  temperature  is  known  (input)  and  the  sky 
exitance  is  calculated  by  an  empirical  equation  as  a function  of  air 
temperature  (input).  Thus,  the  final  equation  becomes 


i  i k x m 4 

• "*"lr  — = a • c • T -CONT.  -SECTOR. -ABSORB  -COS..,. 

2 xx  i xm  j m ilk 

ra  J 

The  total  flux  density  emitted  by  elements  in  layer  i and  absorbed 
by  a particular  mid-element  in  layer  m at  inclination  k can  be  des- 
cribed by 


& 

kxm  _ T 4 

— - — = o-  c -T 

2 x x 

m 

The  total  flux 
inclination  angle  k 


•ABSORB  • 
m 


9 

E CONT.  -SECTOR. 

, j 1 


18 

[ E cos  ! 


density  absorbed  by  a mid-element  in  layer  m at 
is  computed  by  summing  all  sources: 


m x=l  m 

where  x = 1,  2,  3,  4,  5 represents  the  sky.  Layer  1,  Layer  2,  Layer 
3,  and  the  ground,  respectively. 

Nine  equations  for  each  layer  are  constructed.  Each  equation 
represents  the  absorbed  flux  density  for  each  mid-element  inclination. 
For  each  layer  the  appropriate  equation  is  weighted  by  the  frequency 
of  occurrence  of  the  elements  within  the  corresponding  inclination 
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class.  The  nine  equations  are  then  summed  to  represent  the  average 

absorbed  thermal  flux  density  within  the  three  canopy  layers.  In 

2 

addition,  the  flux  density  absorbed  (w/m  ) is  on  a per  unit  area 
2 

basis.  Thus,  the  m“  term  above  must  represent  both  the  top  and 
bottom  surfaces  of  the  leaf.  As  a consequence,  the  factor  of  4 is 
introduced : 

<?  9 


= 41 


k=l  m 


•FKI'QD, 

' km 


— = the  average  absorbed  thermal  flux  density  by  the 

m 

mid-elements  in  layer  m 

FREQD,.^  = the  probability  of  occurrence  of  inclination  k 
for  elements  in  layer  m. 

The  resulting  three  equations  for  each  layer  represent  the  average 
absorbed  thermal  flux  density.  To  complete  the  energy  budget  for 
each  layer  we  must  include:  absorbed  solar  radiation,  convection, 
transpiration,  and  thermal  radiant  emission. 

Solar  Radiation  Absorption 

Several  models  have  been  developed  to  study  the  interactions  of 
solar  radiation  within  vegetation  canopies.  Allen  and  Richardson 
(1968),  Alderfer  and  Cates  (1971),  and  Suits  (1972)  have  adapted  a 
system  of  simultaneous  differential  equations  developed  by  Kubolka 
and  Munk  (1911)  in  various  ways  to  vegetation  canopies.  Suits  (1972) 
developed  a model  which  includes  geometric  effects  and  predicts 
non-Lambor tian  characteristics  of  vegetation  canopies. 
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Another  approach,  developed  by  Oliver  and  Smith  (1974),  is  the 
Solar  Radiation  Vegetation  Canopy  (SRVC)  model.  The  model  is 
stochastic  in  nature  and  predicts  the  diurnal  apparent  directional 
spectral  reflectance  of  a vegetation  canopy.  The  same  geometry 
descriptors  as  described  above  are  utilized  within  the  SRVC  model. 

To  calculate  the  average  absorbed  solar  radiation  within  each 
canopy  layer,  it  is  important  to  include  the  complex  scattering  of 
light  as  a function  of  canopy  geometry.  For  the  purpose  of  calcu- 
lating the  absorbed  solar  radiation,  it  is  believed  that  the  mathe- 
matical framework  of  the  SRVC  model  is  the  most  physically  based  and 
the  most  easily  adapted  to  calculate  the  absorbed  solar  radiation  of 
the  above  models.  A complete  description  of  the  SRVC  model  is 
presented  by  Oliver  and  Smith  (1974). 

The  SRVC  model  has  been  modified  to  estimate  spectral  absorption 
within  vegetation  canopies  and  has  been  specifically  applied  to  a 
cluster  of  four  modeling  lodgepole  pine  trees  (Rimes  and  Smith,  1979). 
The  results  show  that  the  total  global  irradiance  absorbed  by  the 
lodgepole  pine  canopy  system  is  relatively  constant  with  solar  zenith 
angle.  However,  the  proportion  of  total  global  irradiance  absorbed 
by  individual  canopy  layers  varies  as  a function  of  solar  zenith 
angle  (Figure  6).  The  mean  total  solar  flux  absorbed  per  unit  canopy 
element  surface  area  for  any  given  layer  and  solar  zenith  angle  can 
be  readily  calculated  from  the  information  in  Figure  8 as  presented 
by  Rimes  and  Smith  (1979). 


PROPORTION  OF  SOLAR  IRRADIANCE  ABSORBED 


Layer  I 
Layer  2 
Layer  3 
Ground 


SOLAR  ZENITH  ANGLE  ( DCG. ) 


Simulated  proportion  of  global  solar  irrndiance  absorbed 
by  the  lodgepolo  pine  canopy  system  (total).  Layer  1, 
Layer  2,  Layer  3,  and  the  ground,  as  a function  of  solar 
zenith  angle  for  October  14,  1977. 
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Other  Energy  Transfers 

Other  energy  transfers  to  and  from  the  mid-elements  include: 

thermal  exitance,  transpiration,  and  convection. 

The  thermal  exitance  of  all  mid-elements  is  calculated  by  the 

Steffan  Boltzmann  Law: 

4 

M = o • c • T 

i i i 

where: 

M = mean  exitance  of  mid-elements  in  layer  i 

i 

a = Steffan  Boltzmann  constant 

= emissivity  of  mid-elements  in  layer  t (input) 

T ^ = mean  surface  temperature  of  mid-elements  in  layer  i (°K). 
The  ground  thermal  exitance  is  calculated  in  a similar  fashion. 

Gates  (1968)  presented  the  equation  used  for  transpiration,  and 
Lee  and  Gates  (1964)  discussed  the  formulation  in  detail.  The  driving 
force  is  the  difference  between  the  water  vapor  density  within  the 
leaf  and  in  the  free  atmosphere  beyond  the  boundary  layer.  The  water 
vapor  density  within  the  leaf  is  assumed  to  be  at  saturation  at  the 
leaf  temperature.  Controlling  variables  on  transpiration  include  the 
resistance  to  diffusion  offered  by  the  diffusion  pathway,  such  as  the 
stomata  and  the  boundary  layer.  The  equation  for  any  particular  mid- 
leaf is 

so  1 (T  ) - RH*spa(T  ) 

TRANS  = H(T  ) ( — - — ) -(697.8) 

l l K "t-  K 

i a 

where : 

_2 

TRANS^  = transpirational  loss  from  leaf  in  layer  i (w*m  ) 

H ( T ^ ) = latent  heat  of  vaporization  of  water  at  the  leaf 

temperature  T (cal'gm  ^) 


- 
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V 


t 

■t 


i 


spl(T  ) = water  vapor  density  inside  the  leaf  at  saturation 

-3 

at  the  leaf  temperature  (gm*cm  ) 


RU 


= relative  humidity  of  air  (input) 


spa(T  ) = water  vapor  density  at  saturation  of  the  free  air 

cl 

beyond  the  boundary  layer  of  the  leaf  at  the  air 

-3 

temperature  T (gn*cm  ) 

cl 

R = internal  leaf  resistance  to  water  vapor  diffusion 

i 

in  min* cm  ^ (input) 


R = resistance  of  the  boundary  layer  to  water  vapor 

cl 

diffusion  (min*cm  ^). 

H(T  ),  spl(T  ),  and  spa(T  ) were  calculated  using  physically  based 

l l cl 

formulas.  Values  of  R for  lodgepole  pine  needles  were  estimated 

cl 

from  a mass  transfer  determination  of  Landsberg  and  Ludlow  (1969)  who 
used  Sitka  spruce  shoots.  The  formula  is 


R.  = (0.04  + 1.27  (p-°'5))/60 

d 


where: 

P = wind  speed  in  cmvsec  ^ (input). 

The  constant  R^  value  used  for  the  lodgepole  pines  in  this  study 
was  0.66  min/cm.  Gates  (1966)  and  Miller  and  Gates  (1967)  reported 
R^  values  of  0.72,  0.33,  and  0.50  min/cra  for  Picea  mariana,  Pinus 
res inosa,  and  Pinus  strobus,  respectively.  Jarvis  et  al . (1976)  and 
Tenhunen  and  Gates  (1975)  presented  recent  investigations  of  the 
stomatal  opening  and  closing  as  influenced  by  environmental  factors 
and  concluded  that  the  complex  control  of  the  stomata  lias  not  yet 
been  described  adequately. 


i 
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Since  the  mid-elements  represent  both  leaves  and  branches  when 
dealing  with  woody  plants,  one  can  assume  that  branches  do  not  trans- 
pire and  the  transpiration  equation  can  be  weighted  according  to  the 
branch  area  index  and  the  leaf  area  index. 

The  following  convection  equation  was  utilized.  Tibbals  et  al . 
(1964)  conducted  quantitative  measurements  on  silver  castings  of  blue 
spruce  and  white  fir  branches  in  a controlled  radiation  and  windtunnel 
chamber.  The  authors  report  convective  coefficients  for  free  convec- 
tion in  both  species.  However,  Gates  (1968)  notes  that  rarely  in 
nature  is  there  any  air  movement  less  than  S.S  cm "sec  ^ (0.2  m.p.h.). 
As  a consequence,  an  equation  describing  forced  convection  can  be 
used  to  approximate  all  convectional  exchanges.  Tibbals  et  al^.  (1964) 
found  that  both  longitudinal  and  horizontal  wind  flows  gave  equal 
coefficients  for  spruce. 

For  P>30.0  h = (0.95P0,97) • (.698) 

c 

For  jj < 30 . 0 h = (20.4  + 0. 2g° ' 97)  • ( . 698) 

c 

where: 

p = wind  velocity  in  cm*sec  ^ (input) 

h = convectional  coefficient  in  wm  7-°C 
c 

The  convectional  exchange  of  a mid-element  is  calculated  as 


QFC  = hc 


(T  - T ) 
s a 


where: 


h 

c 


T 

s 

T 

a 


power  per  unit  area  of 
convection  coefficient 
surface  temperature  of 
air  temperature  of  the 
layer  in  °C  (input). 


mid-element  loss  or  gain  (wm  “) 

. -2  0 -1 
in  w*m  • C 

mid-element  in  °C 

free  air  beyond  the  boundary 
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Kind  speed  is  highly  variable  from  point  to  point  within  the  canopy 
(Bergen,  1974).  As  a consequence,  the  mean  measured  wind  speed 
values  were  utilized  for  all  three  canopy  layers. 

The  sky  thermal  exitance  was  calculated  by  an  empirical  equation 
dependent  only  on  air  temperature  near  the  ground  surface,  and  clear 
sky  conditions  were  assumed.  Hudson  (1969)  presented  several 
references  which  estimate  sky  thermal  exitance. 

It  is  imporcant  to  note  that  a multitude  of  convectional,  trans- 
pirational,  sky  thermal  exitance,  and  solar  absorption  formulations 
exist  that  may  be  more  suitable  for  specific  modeling  objectives. 

For  this  reason  the  model  has  been  structured  so  that  Afferent  for- 
mulations of  the  above  can  be  easily  incorporated  within  the  model. 


Model  Solution 

The  total  energy  budget  equations  for  each  canopy  layer  can  now 
be  formed.  The  result  is  a system  of  three  nonlinear  equations  and 
three  unknowns  being  the  surface  temperature  of  the  mid-elements  in 
each  layer  which  represent  the  respective  average  temperature  of 
each  layer. 

[Layer  1 energy  budget  equation] 


-> 


F = Layer  2 energy  budget  equation  = 0 
Layer  3 energy  budget  equation 
To  solve  this  system  of  equations  the  model  calls  the  ZSYSTM 
algorithm  which  exists  in  the  International  Mathematical  and  Statis- 
tical Library  ( 1° 7 7 ) . ZSYSTM  is  a numerical  routine  which  uses 
Brown's  method  (1969,  1971)  for  solving  N simultaneous  nonlinear 
equations  in  N unknowns.  The  method  is  at  least  quadratic  convergent 
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2 

and  requires  only  N /2  + 3N/2  function  evaluations  per  iterative  step 
2 

as  compared  with  N + N evaluations  for  Newton's  method. 

The  roots  of  the  system  predict  the  average  temperature  of  the 
layers  and  are  used  to  calculate  the  following  thermal  predictions. 


Thermal  Predictions 

The  model  predicts  the  thermal  radiance,  effective  radiant 
temperature  (ERT) , and  equivalent  exitance  in  the  9 viewing  inclina- 
tion bands  at  10°  intervals  above  the  canopy.  The  contribution  of 
each  layer  and  the  ground  to  the  nine  sensing  positions  are  calcu- 
lated as  follows.  The  thermal  radiance  in  the  band  directions  j are 

-1  4 

L.  = it  • [ PHIT . , *e  •a*X1 
J jl  1 1 


-f. 

PGAP . . 

•PHIT.  • 

,,  4 

c. -a*X„ 

jl 

J2 

2 2 

+ 

PGAP . . 

•PGAP . „ • 

PHIT  -e. -a-X 

Jl 

j2 

J 3 3 

+ 

PGAP . , 

•PGAP  • 

PGAP .•£. -a-X 

Jl 

J 2 

J 3 4 

where: 

each  row  represents  the  thermal  radiance  contribution  to 
the  sensor  by  Layer  1,  Layer  2,  Layer  3,  and  Layer  4 (ground), 
respectively . 

Lj  = thermal  radiance  of  the  sensor  at  viewing  angle  j 

, -2  -1. 
i (wm  *sr  ) 

X = average  surface  temperature  of  elements  in  laver  m (°K): 
ra  1 

m = 1 , 2 , 3 , 4 . 

The  thermal  radiance  (Lj)  can  be  converted  to  the  equivalent 
exitance  (M.)  by 


M.  = L *tt 

j j 
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and  the  effective  radiant  temperature  (EKT\°K)  in  band  direction  j 
can  be  calculated  as 


where  a is  the  Stefan-Boltzmann  constant. 

The  model  also  predicts  the  response  of  a thermal  sensor  looking 
horizontally  from  the  ground  into  any  of  the  three  layers.  When 
looking  horizontally  at  a canopy  the  probability  of  gap  is  0.0  accord- 
ing to  the  assumptions  in  the  model.  Thus,  for  a relatively  narrow 
field  of  view  the  ERT  of  any  given  layer,  when  looking  horizontally, 
is  calculated  simply  by  using  the  Stefan-Boltzmann  equation  with  the 
appropriate  emissivity  factor  and  average  layer  temperature. 

In  addition,  for  each  simulation,  the  average  predicted  temper- 
ature of  each  layer,  ground  thermal  exitance,  sky  thermal  exitance, 
absorbed  solar  flux  density  of  each  layer,  thermal  exitance  of  each 
layer,  absorbed  thermal  flux  density  of  each  layer,  convectional 
exchange,  transpirational  exchange,  the  geometric  coefficients 
(including  CONT,  COS,  SECTOR,  LAI,  S,  FREQD)  and  all  input  parameters 
are  displayed. 

Field  Measurements 

A unique  thermal  and  environmental  data  base  for  a lodgepole 

pine  canopy  at  Leadvilla,  Colorado  was  collected  during  1977.  Four 

clustered,  lodgepole  pine  trees  were  chosen  for  intensive  study  as 

shown  in  Figure  7.  These  modeling  trees  had  the  following  mean 

statistics:  6.0  m height,  30  year  age,  13.2  cm  DBH,  and  a surrounding 
2 

stand  of  102  m /hectare  basal  area.  The  S parameter,  foliage  area 

4 I 
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indices,  and  foliage  angle  frequency  distributions  of  the  modeling 
canopy  were  measured  as  reported  by  Kimes,  Smith,  and  Berry  (1979), 
and  Kimes  and  Smith  (1979). 

Personnel  from  the  Army  Corps  of  Engineers  Waterways  Experiment 
Station,  Environmental  Laboratory  (WES/EL)  at  Vicksburg,  Mississippi 
developed  a system  for  automated  collecting,  processing  and  displaying 
environmental  baseline  data  as  described  by  West  and  Floyd  (1976). 

The  system  was  utilized  to  monitor  environmental  conditions  at  the 
study  site  for  the  months  of  July,  September,  and  October,  1977.  All 
sensor  measurements  were  recorded  once  every  hour  continuously  for  the 
duration  of  the  study.  The  measurements  taken  included  air  temper- 
ature, global  solar  irradiance,  wind  speed,  wind  direction,  rainfall, 
soil  temperature,  and  vegetation  surface  temperature.  The  specific 
make  and  calibration  procedures  of  the  above  instrumentation  are 
described  by  West  and  Floyd  (1976). 

Figure  7 is  an  oblique  photograph  of  the  study  area  showing  the 
sensor  positions.  At  station  Ml  in  the  meadow  opening,  air  tempera- 
ture, wind  speed,  wind  direction,  global  solar  irradiance,  and  precip- 
itation were  measured.  Air  temperature  was  measured  at  a height  of 
1 m.  All  of  the  above  measurements  were  also  taken  at  station  M2 
within  the  tree  area.  In  addition,  the  air  temperature  within  the 
center  of  the  four  modeling  trees  (M3)  was  recorded. 

The  surface  temperature  of  the  base  of  needles  at  three  branch 
tips  of  the  modeling  trees  were  measured  by  contact  thermistors;  one 
thermister  was  located  1 m above  the  ground  in  Layer  3,  the  second 
thermistor  was  located  2 m above  the  ground  in  Layer  2,  and  the  final 
thermister  was  located  in  the  top  of  the  modeling  trees. 


Oblique  photograph  of  modeling  trees,  the  meteorological 
measurement  stations  (Ml,  M2,  M3)  and  the  4 stake  positions 
(SI,  S2,  S3,  S4) . 
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Several  other  supporting  data  were  recorded.  During  July  15-16, 
1977,  the  Wahl  Digital  Heat  Spy  -DSH-  14  Thermal  Radiometer  with  a 
band  pass  of  4.8  to  20.0  Pm  and  a 3.5°  field  of  view  was  used  to 
measure  the  average  horizontal  effective  radiant  temperature  (ERT)  of 
Layer  2 of  the  canopy  at  four  stake  positions  (Figure  7).  Simulta- 
neously, black  and  white  Polaroid  photographs  derived  from  the  AGA- 
Thernovision  were  taken.  The  system  operates  in  the  2.0  to  5.6  pm 
region  and  scans  a particular  scene.  In  addition,  during  August  14-15, 
the  Wahl  Heat  Spy  was  used  to  measure  needle  temperature  at  five 
branch  tips  in  the  modeling  trees. 

Data  Reduction 

An  initial  analysis  of  the  Wahl  Heat  Spy  and  contact  thermister 
data  was  completed  to  decide  which  data  form  was  most  suitable  for 
testing  the  TCSM's  accuracy  of  prediction. 

The  Heat  Spy  and  contact  thermister  were  taken  of  the  needle- 
branch  tip  complex  which  has  a relatively  large  concentration  of  mass 
(Figure  8).  Due  to  the  small  diameters  of  the  needles  it  was  not 
feasible  to  securely  fasten  the  contact  thermistors  to  them.  As  a 
consequence,  the  thermisters  were  placed  between  the  bases  of  the 
needles  at  the  branch  tips.  The  bases  of  the  needles  at  the  branch 
tips  were  very  concentrated  and  thus,  completely  surrounded  the 
contact  thermisters  (Figure  8). 

The  Heat  Spy  ERT's  were  correlated  with  the  temperatures  of  the 
contact  thermisters  on  many  different  natural  and  man-made  surfaces 
during  the  month  of  July  at  Vicksburg,  Mississippi  by  WES/EL  person- 
nel. The  results  showed  that  the  two  instruments  corresponded  within 
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FOV  of  heat  spy  on  moss  of  needle  bases  at  branch  tip. 


BRANCH  TIP 


r 


Figure  8.  Diagram  of  a branch 
Heat  Spy  Radiometer 
thermis  ter. 


4 


tip  showing  the  target  of  the  Wahl 
and  the  placement  of  the  contact 


1°C  under  homogeneous  solar  heating  loads.  Thus,  it  was  believed  that 
the  two  measurements  correlated  well  for  the  same  target. 

Two  branch  tips  on  the  modeling  trees  had  contact  thermisters 
placed  at  the  base  of  the  needles  next  to  the  branch  tip  and  were  also 
simultaneously  measured  with  the  Heat  Spy.  In  general,  the  contact 
thermisters  recorded  higher  temperature  values,  of  2°C  or  greater, 
relative  to  the  Heat  Spy  ERT's.  One  possible  explanation  of  this 
deviation  in  measurements  is  as  follows. 

The  contact  thermisters  were  buried  in  a relatively  large  mass 
of  compressed  needles  which  occurs  near  the  branch  tip.  Under  this 
situation,  relatively  little  convectional  or  transpirational  exchange 
of  these  compressed  portions  of  the  needles  can  occur,  which  may 
account  for  the  discrepancy  seen  above.  In  addition,  some  conduction 
from  the  warmer  branches  may  be  operating  as  will  be  discussed  later. 

In  contrast,  the  field  of  view  of  the  Heat  Spy  was  on  the  branch 
tip  (Figure  S)  which  was  exposed  to  the  air.  It  is  believed  that  this 
portion  of  the  branch  tip  undergoes  limited  convectional  and  trans- 
pirational exchange.  As  a result,  the  Heat  Spy  ERT's  were  lower  than 
those  of  the  corresponding  contact  thermisters. 

In  this  modeling  study  the  criteria  used  for  validation  is  the 
accuracy  of  prediction  of  the  mean  canopy  element  surface  temperatures 
and/or  the  mean  horizontal  ERT  of  a canopy  layer.  Gary  (1976)  has 
shown  that  the  needle  area  index  for  a particular  lodgepole  pine  tree 
accounts  for  approximately  87%  of  the  total  element  surface  area  of 
the  tree.  The  large  majority  of  the  needles  extend  into  free  space 
and  undergo  relatively  high  convectional  and  transpirational  exchange 
due  to  the  small  needle  diameters.  As  a consequence,  the  majority  of 
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the  canopy's  surface  area  should  closely  approximate  air  temperature 
(Gates,  1975). 

The  above  measurements  (the  contact  thermister  temperatures  and 
Heat  Spy  ERT's)  of  the  branch  tip  do  not  reflect  the  surface  temper- 
atures of  the  needles  extending  into  free  space,  and  bias  the  measured 
canopy  temperature  to  be  too  high.  Therefore,  these  measurements  are 
generally  above  air  temperature  and  do  not  reflect  the  mean  canopy 
element  surface  temperatures.  Initial  comparisons  between  the  simu- 
lated mean  element  surface  temperatures  for  the  three  layers  for 
October  14-15  and  the  three  contact  thermister  measurements  in  Layers 
1,  2,  and  3 show  a very  poor  accuracy  of  prediction. 

It  is  believed  that  the  horizontal  ERT's  as  measured  by  the  Wahl 
Heat  Spy  from  the  four  stake  positions  during  July  15-16  were  the 
least  biased  of  all  validating  measurements,  since  the  field  of  view 
incorporated  a cross  section  of  canopy  element  types.  However,  during 
the  August  14-15  date  only  the  ERT  of  five  branch  tips  were  taken. 

Simulations 

The  thermal  behavior  of  the  modeling  canopy  was  simulated  for 
two  complete  diurnal  cycles.  For  the  July  15-16,  1977  simulation, 
the  mean  of  the  Heat  Spy  ERT's,  as  measured  from  the  4 stake  positions, 
was  used  to  test  the  accuracy  of  prediction.  Two  simulations  using 
the  air  temperatures  within  the  canopy  (M3)  and  in  the  meadow  opening 
(Ml)  were  performed  during  this  date. 

One  simulation  was  performed  for  October  14-15,  1977.  The  Heat 
Spy  measurements  of  branch  tip  needles  were  utilized  to  test  the 
accuracy  of  prediction. 

i 
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Sensitivity  analyses  were  performed  on  the  following  model 
parameters:  average  canopy  element  emissivity,  average  canopy  element 
internal  resistance  to  water  vapor  diffusion,  and  canopy  geometry. 

Res ults  an d_ Discuss  ion 

The  simulated  horizontal  ERT's  for  the  three  layers  for  July 
15-16,  using  the  air  temperature  probe  in  the  middle  of  Layer  2 of 
the  modeling  canopy  (M3),  are  presented  in  Figure  9 along  with  the 
mean  horizontal  Heat  Spy  ERT's  from  the  A stake  positions.  The 
corresponding  measured  solar  irradiance  and  air  temperature  are 
presented  in  Figure  10.  Wind  speed  was  0.0  m/s  for  all  measurement 
periods.  The  minimum  recorded  wind  speed  possible  was  10  cm/s. 

Gates  (1968)  stated  that  rarely  is  wind  speed  in  natural  environments 
below  S cm/s.  As  a consequence,  for  periods  when  wind  speed  was 
recorded  to  be  0.0,  a minimum  value  of  10  cm/s  was  utilized.  It 
should  be  noted  that  the  simulated  data  were  derived  from  the  meteor- 
ological data  which  were  recorded  at  hourly  intervals  and  the  Heat 
Spy  measurements  were  not  necessarily  synchronous  in  time.  Conse- 
quently, one  must  compare  the  general  trends  of  the  simulated  data 
with  the  Heat  Spy  measurements.  In  fact,  the  erratic  nature  of  the 
solar  irradiance  (Figure  10)  suggests  that  between  hourly  intervals, 
the  true  solar  irradiance  function  could  vary  widely.  This  fact 
could  explain  some  of  the  deviations  during  the  day  shown  in  Figure  9. 
During  the  night  the  simulated  values  deviated  from  the  measured 
temperature  by  less  than  1.5°C. 

Selected  output  for  0930  and  0330  (Standard  Time)  simulations 
is  presented  in  Table  1.  During  the  day  the  average  layer  temperature 
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» Layer  I Simulated 
O Layer  2 Simulated 
A Layer  3 Simulated 
Q Heat  Spy  Measurement 


Figure  9. 


STANDARD  TIME 


Simulated  versus  measured  lodgepole  pine  canopy  horizontal 
ERT's  for  July  15-16,  1977.  Measured  ERT's  are  the  mean 
of  4 horizontal  ERT's  of  the  middle  layer  as  measured  by 
the  Wahl  Heat  Spy.  Air  temperature  was  recorded  within 
center  of  canopy  (M3  site) . 
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Table  1.  Selected  output  for  0930  and  0330  (Standard  Time) 
July  15-16,  1977. 


Time  = 0930 

Average  element  temperatures  (Layers  1-3)  = 20.4,  16.6,  16.3°C 

The  thermal  exitance  and  ERT  above  the  canopy  for  the  various  viewing 
angles  are: 


Inc lination  (degrees) 


Exitance  (v/m  ) 


ERT  (°C) 


5 

419 

20.3 

15 

413 

19.2 

25 

410 

13.6 

35 

408 

18.4 

45 

408 

18.2 

55 

407 

18.1 

65 

407 

18.1 

75 

407 

18.1 

85 

407 

18.1 

Time  = 0330 

Average  element  temperatures  (Layers  1-3)  = -1.0,  -0.0,  0.5°C 

The  thermal  exitance  and  ERT  above  the  canopy  for  the  various  viewing 
angles  are: 


Inclination  (decrees) 

9 

Exitance  (w/m  ) 

ERT  (°C) 

5 

310 

-1.0 

, 15 

311 

-0.6 

25 

313 

-0.4 

35 

314 

-0.2 

45 

314 

-0.1 

55 

314 

-0.1 

65 

315 

-0.1 

75 

315 

-0.1 

85 

315 

-0.1 

\ 

( 
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decreases  as  one  proceeds  from  Layer  1 to  Layer  3 due  to  solar  heating 
and  canopy  geometry  interactions.  During  the  nighttime,  however,  the 
layers  cool  differentially  due  to  the  relatively  low  thermal  exitance 
of  a clear  sky  and  the  relatively  high  surface  temperature  of  the 
ground.  As  a consequence.  Layer  3 has  the  highest  temperature.  The 
three  contact  thermisters  in  Layers  1,  2,  and  3 support  this  relation- 
ship. These  trends  are  discussed  by  Geiger  (1961). 

The  ERT  of  a sensor  above  the  canopy  as  a function  of  view  angle 
is  dependent  on  the  above  layer  temperature  differentials  and  the 
canopy  geometry.  At  the  lowest  sensor  inclination  angles,  the  ERT 
strongly  relfects  the  temperature  of  Layer  1 (Table  1).  As  the  sensor 
view  inclination  angle  increases,  the  second,  third,  and  ground  layer 
temperatures  more  strongly  influence  the  sensor  ERT.  It  is  believed 
that  particular  canopy  geometries  (leaf  angle  distribution,  leaf  area 
index,  and  leaf  spatial  distributions)  can  have  very  significant 
effects  on  the  ERT  of  the  sensor  at  varying  view  angles. 

Sensitivity  analyses  were  run  at  two  different  times  (0930  and 
0330)  during  July  15-16  to  characterize  two  extremes  of  environmental 
conditions.  The  fixed  parameters  for  Layers  1,  2,  and  3 used  in 
conducting  the  sensitivity  analysis  are  as  follows.  The  emissivity 
(including  the  ground),  wind  velocity,  and  leaf  resistance  to  water 
vapor  diffusion  were  1.0,  10.0  cm/s,  and  0.66  min/cm,  respectively, 
for  both  0930  and  0330  time  periods.  The  solar  irradiance,  air  temp- 
erature, ground  temperature  and  relative  humidity  for  the  0930  time 
2 

were  855  w/m  , 14.6°C,  11.7°C,  and  0.20,  respectively,  and  for  the 

2 

0330  time  were  0.0  w/m  , 0.4°C,  5.0°C,  and  0.85,  respectively. 


k 


i 
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The  results  of  the  sensitivity  analysis  of  emissivity  for  the 
day  and  night  environmental  conditions  show  that  within  a reasonable 
emissivity  range  for  natural  vegetation  (0.96-1.00)  the  change  in 
average  element  temperature  is  on  the  order  of  0.6°C  for  all  layers. 


\ change  of  the  internal  resistance  to  waver  vapor  diffusion  parameter 

-1 


(R^)  within  a range  of  0. 3-1.2  min*cm  has  an  equally  small  effect 
on  the  average  element  temperature  of  the  three  layers  for  the  day 
environmental  conditions.  However,  at  lower  values  the  parameter  is 
very  sensitive  for  day  environmental  conditions.  For  the  night  condi 


tions,  the  average  element  temperature  of  the  three  layers  changed 
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less  than  0.3°C  for  the  range  of  0.05-1.2  min*cm 


In  this  study  the  R^  was  held  constant  at  0.66  min*cm  . Running 


(1978)  has  found  that  is  variable  for  lodgepole  pine  needles  in 


full  sunlight.  The  minimum  daily  R^  may  vary  from  0.11  to  0.50 


min'cm  depending  on  the  pre-dawn  leaf  water  potential.  The  daily 
variation  of  R^  is  largely  dependent  on  the  humidity  and  may  vary  as 


much  as  four  times  that  of  the  minimum  daily  R^.  Thus,  the  assumption 


of  a constant  R^  is  very  erroneous.  As  the  season  progresses,  the 
pre-dawn  leaf  water  potential  drops  as  a result  of  decreasing  avail- 


able soil  water,  and  during  the  October  simulation  a constant  R^  value 


-1 


of  0.30  min*cm  or  above  would  be  a relatively  good  value  for  both 
day  and  night  conditions  in  light  of  the  fact  that  the  insensitivity 


of  R^  on  canopy  temperatures  was  above  this  value.  However,  during 


the  July  simulation  for  full  sunlight  on  the  needles,  an  R^  value  on 


-1  . 


the  order  of  0.15  min-cm  is  more  appropriate.  Using  the  day  environ- 


mental conditions,  the  average  element  temperature  in  Layer  1 decreases 


-1 


3 C for  an  change  from  0.66  to  0.15  min*cm  . This  tendency  would 


‘ ,*  vjr  • k 
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cause  the  temperature  of  Layer  1 to  be  much  closer  to  air  temperature 
for  the  July  simulation  (Figures  9 and  10).  R^  of  needles  tend  to  be 
maximal  at  radiation  value,  less  than  10%  of  full  sunlight,  and  sto- 
mata are  generally  closed  in  the  dark  (Hinckley,  Lassoie,  and  Running, 
1978).  Thus,  in  the  relatively  shaded  portions  of  the  canopy  (Layers 
2 and  3)  one  would  expect  the  average  R to  be  high,  and  Layer  1, 
which  intercepts  a large  proportion  of  the  solar  irradiance,  would 
have  a minimal  R value. 

i 

Tan,  Black,  and  Nnyamah  (1978)  have  noted  that  in  well  ventilated 
coniferous  canopies,  leaf  temperatures  are  relatively  similar  to  air 
temperature.  Figures  9 and  10  indicate  the  close  correspondence  be- 
tween air  temperature  and  the  simulated  roan  cr.  . pv  element  tempera- 
tures for  Layers  2 and  3.  However,  Layer  1 significantly  deviates 
from  air  temperature  during  times  of  high  solar  irradiance.  If  more 
appropriate  values  of  for  each  layer  would  have  been  used,  as 
suggested  above,  the  simulated  average  layer  temperatures  would 
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approximate  air  temperature. 

Table  2 presents  the  average  element  temperature,  sensor  ERT, 
and  probability  of  gap  for  two  canopy  geometries,  and  for  both  day 
and  night  environmental  conditions.  The  two  geometries  used  were  the 
normal  foliage  angle  distribution,  which  was  measured  and  utilized  in 
all  other  analyses  above,  and  the  theoretical  erectophile  foliage 
angle  distribution.  Both  distributions  are  presented  by  Rimes,  Smith, 
and  Berry  (1978).  It  should  be  noted  that  although  canopy  geometry 
clearly  affects  the  manner  in  which  solar  radiation  is  absorbed  by  a 
canopy  system,  and  thus  the  layer  temperature,  the  solar  flux  absorbed 
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by  each  layer  was  hold  constant  for  this  sensitivity  analysis 
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Table  2.  Sensitivity  analysis  for  the  effect  of  two  canopy  geometries 
(normal  and  erectophile)  versus  average  element  temperature 
(A)  and  effective  radiant  temperature  (ERT) (B)  above  the 
canopy  as  a function  of  view  angle  for  both  day  and  night 
environmental  conditions.  In  addition,  the  probability  of 
gap  (PGAP)  for  the  nine  i elination  intervals  and  each 
layer  are  compared  between  the  two  canopy  geometries  (C) . 


(A)  Average  Element  Temperature  (°C) 


Layer 

Day 

Night 

Normal 

Erectophile 

Normal 

Erectophile 

1 

20.4 

20.8 

-1.0 

-0.8 

2 

16.6 

16.6 

0.0 

0.0 

3 

16.3 

16.4 

0.5 

0.4 

(B)  ERT 

(°C) 

Inclination 

Day 

Night 

angle 

Normal 

Erectophile 

Normal 

Erectophile 

5 

20.3 

20.8 

-1.0 

-0.8 

15 

19.2 

20.3 

-0.6 

-0.7 

25 

18.6 

19.5 

-0.4 

-0.5 

35 

18.4 

18.7 

-0.2 

-0.2 

45 

18.1 

18.0 

-0.1 

0.2 

55 

18.1 

17.4 

-0.1 

0.6 

65 

18.1 

16.8 

-0.1 

1.0 

75 

1S.1 

16.3 

-0.1 

1.4 

85 

1S.1 

16.0 

-0.1 

1.6 

(C)  PGAP 


Inclination  Layer  1 Layer  2 Layer  3 


angle 

Normal 

Erectophile 

No  rma 1 

Erectophile 

Normal 

Erectophile 

5 

.03 

.00 

.00 

.00 

.13 

.00 

15 

.30 

.11 

.10 

.02 

.47 

.23 

25 

.41 

.30 

.19 

.10 

.57 

.44 

35 

.46 

.44 

.23 

.21 

.61 

.57 

45 

.48 

.54 

.25 

.31 

.62 

.66 

55 

.49 

.62 

.26 

.40 

.63 

. 72 

65 

.50 

.68 

.27 

.48 

.63 

.77 

75 

.50 

.72 

.27 

.54 

.63 

.80 

85 

.50 

.74 

.27 

.57 

.63 

.82 

• • s. 
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As  shown  in  Table  2-A,  the  effect  of  canopy  geometry  on  thermal 
radiation  transfers  is  minimal  and  as  a consequence  the  average  ele- 
ment temperature  changes  very  little  between  the  two  canopy  geometries. 
However,  the  canopy  geometry  clearly  has  an  effect  on  the  contribution 
of  thermal  radiation  from  each  layer  to  the  sensor  above  the  canopy 
as  a function  of  view  angle  (Table  2-B) . These  trends  can  be  explained 
by  the  different  probabilities  of  gap  at  the  various  view  angles  for 
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each  layer  as  seen  in  Table  2-C. 

The  simulated  results  in  Figure  9 were  derived  from  the  air 
temperature  probe  in  the  center  of  the  four  modeling  trees  (M3) , and 
the  correspondence  between  simulated  and  measured  data  was  relatively 
good.  However,  it  was  noted  that  individual  air  temperature  probes 
were  not  well  correlated.  For  example,  the  air  temperatures  as 
measured  by  the  air  temperature  probes  at  the  Ml,  M2,  and  M3  sites 
were  compared  for  the  July  14  and  16  diurnal  cycles.  During  the  day- 
light hours  the  individual  probes  were  highly  erratic  and  uncorrelated, 
and  at  night  the  probe  in  the  meadow  opening  (Ml)  was  consistently 
lower  than  the  probes  within  the  canopy  (M2,  M3)  by  as  much  as  2°C. 

As  a consequence,  the  simulated  results  for  July  15-16,  using  the  air 
temperature  in  the  meadow  opening  (Ml),  were  relatively  poor,  espe- 
cially at  night,  in  accuracy  of  prediction  as  compared  to  Figure  9. 

The  above  2°C  temperature  differential  at  night  can  explain  these 
relatively  poor  results. 

Bergen  (1971,  1974)  showed  that  air  temperature  differentials 
within  lodgepole  pine  canopies  can  be  as  much  as  4-5°C  in  the  vertical 
profile  for  clear,  sunny  days.  In  addition,  six  simultaneous  air 
temperature  measurements  at  various  horizontal  points  within  the 
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canopy  shoved  air  temperature  differences  as  large  as  2°C.  Thus,  air 
temperature  measurements  at  a single  point  (1  m above  the  ground)  as 
utilized  within  this  study  will  introduce  error.  In  addition,  on 
cloudless  nights,  ground  cooling  by  net  radiation  loss  often  occurs, 
and  an  air  temperature  inversion  near  the  ground  occurs.  Geiger 
(1961)  has  shown  that  even  minimal  wind  speed  of  10  to  100  cm/s  can 
disrupt  temperature  stratification  during  the  night  in  which  2°C 
differentials  can  commonly  occur.  These  fluctuations  can  introduce 
error . 

In  addition,  Bergen  (1971)  presented  windspeed  variation  for  a 
typical  clear  day  and  at  a typical  station  within  a 10  m tall  lodge- 
pole  pine  stand.  Average  differentials  between  simultaneous  profiles 
were  about  10%  below  the  live  canopy  and  20%  in  the  live  canopy  region. 
A subcanopy  maximum  occurs  near  3 m height  and  a region  of  minimum 
wind  speed  occurs  near  6 m height.  Thus,  some  error  is  introduced  by 
assuming  a homogeneous  vertical  wind  profile.  Jarvis  et  al.  (1976) 
discussed  other  studies  of  wind  speed  profiles  in  conifer  canopies. 

The  data  derived  from  the  AGA  Thermovision  on  July  15-16  were  not 
reduced  to  absolute  temperatures  due  to  several  technical  difficulties. 
However,  the  black  and  white  Polaroid  photographs  derived  from  the 
system  were  utilized  to  document  relative  trends  which  occurred  within 
the  lodgepole  pine  canopy.  The  photographs  show  several  interesting 
trends.  Figure  11-A  shows  the  canopy  at  0700.  Some  of  the  needles 
and  small  branches  tend  to  heat  up  due  to  solar  heating;  however,  this 
phenomenon  is  very  heterogeneous  in  nature  due  to  the  heterogeneous 
distribution  of  sun-flecks  within  the  canopy  as  discussed  by  Ross 
(1976).  This  phenomenon  is  supported  by  the  erratic  contact  thermister 


measurements  of  the  three  branch  tips  in  Layers  1,  2,  and  3,  respec- 
tively, during  the  day.  In  addition,  the  bole  and  larger  branches 
tend  to  be  relatively  cool  due  to  the.  relatively  large  mass  (Figure 
11-A) ; and  thus,  the  steady-state  assumption  of  energy  exchange  may 
not  be  appropriate  for  these  elements.  This  heterogeneity  is  impor- 
tant, especially  in  the  design  of  the  field  measurements  used  to 
validate  the  simulated  results.  The  AGA  photographs  show  at  0900 
that  the  boles  and  larger  branches  are  still  relatively  cooler  than 
the  other  canopy  elements.  During  the  morning  the  above  trends  are 
supported  by  Heat  Spy  measurements.  The  Heat  Spy  measurements  of  4 
branch,  7 bole,  and  5 needle-branch  tip  samples  for  0800  and  1000 
show  mean  values  of  the  branch  and  bole  FRT  to  be  1-4  °C  cooler  than 
the  mean  needle-branch  tip  FRT. 

However,  by  1100  the  branches  (Figure  11-B)  generally  tended  to 
be  warmer  than  the  other  canopy  elements.  Figure  11-C  shows  a close- 
up  of  a group  of  branches  at  1300.  The  branches  generally  have  a 
higher  total  solar  absorption  coefficient  (Rimes  and  Smith,  1979), 
and  they  do  not  possess  the  high  degree  of  conventional  exchange  and 
transpirational  exchange  that  the  needles  experience  due  to  their 
larger  mass  and  physiology.  During  the  afternoon  the  Heat  Spy  measure- 
ments of  4 branch,  7 bole,  and  5 needle-branch  tip  samples  for  1200 
and  1400  show  mean  values  of  the  branch  FRT  to  be  6 to  8°C  warmer  than 
the  needle-branch  tip  FRT.  And  the  mean  bole  ERT  was  2°C  cooler  and 
1°C  warmer  than  the  needle-branch  tip  ERT,  respectively,  for  the  two 
times . 

Figure  11-D  shows  that  in  the  early  morning  hours  the  boles,  which 
have  a high  heat  capacity,  are  still  relatively  warm.  The  ERT  of  the 
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mean  branch  and  needle-branch  tip  samples  were  equal  and  the  mean 
bole  ERT  was  3°C  warmer. 

The  majority  of  photographs  derived  from  the  AGA  Thermovision 
did  not  include  the  uppermost  crowns.  Thus,  it  is  difficult  to  com- 
pare trends  of  mean  layer  temperature  as  a function  of  height  within 
the  canopy  for  both  simulated  and  measured  data.  The  photographs 
largely  emphasize  Layers  2 and  3 of  the  canopy.  As  can  be  seen  in 
Figure  11,  the  simulated  results  for  Layers  2 and  3 are  very  close  in 
absolute  temperature,  and  the  photographs  (Figure  11-A,B)  do  not 
demonstrate  any  clear  trends  of  mean  element  temperature  as  a function 
of  height. 

At  night  the  simulated  results  (Figure  9)  suggest  that  under 
clear  sky  conditions  the  mean  temperature  of  Layer  1 will  be  approx- 
imately 2°C  cooler  than  the  other  layers  due  to  a high  net  thermal 
radiant  loss  to  the  sky.  However,  several  photographs  taken  at  night 
tend  to  show  a portion  of  Layer  1.  A trend  of  cooler  canopy  elements 
with  increasing  canopy  height  can  be  seen  in  Figure  11-D. 

The  1 1-A  photograph  suggests  another  source  of  error.  The  TCSM 
assumes  canopy  layers  of  infinite  extent;  however,  in  the  cluster  of 
four  modeling  trees,  solar  heating  of  the  edges  of  the  canopy  system 
does  occur. 

All  simulated  results  for  July  15-16  were  validated  against  the 
mean  Heat  Spy  ERT  measurements  of  the  4 stake  positions.  It  was 
believed  that  these  measurements  were  the  best  ground  truth  available 
in  regards  to  the  modeling  criteria.  However,  these  measurements 
were  not  taken  for  the  October  14-15  date,  and  the  mean  ERT  of  5 
needle-branch  tips  (Figure  8)  was  utilized  for  validation. 

4 I 
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The  simulated  mean  element  surface  temperatures  for  October  14- 
15  are  presented  in  Figure  12  along  with  the  mean  Heat  Spy  ERT's  from 
5 branch  tips  occurring  in  Layers  1 and  2.  The  global  ir radiance, 
air  temperature,  and  wind  speed  measured  at  the  Ml  site  are  presented 
in  Figures  13  and  14,  3‘he  accuracy  of  prediction  is  relatively  good 
except  at  night  when  a deviation  of  approximately  3°C  consistently 
occurs.  A possible  explanation  of  this  deviation  is  as  follows. 

Only  the  air  temperature  probe  in  the  meadow  opening  (Ml  site) 
was  available  during  the  October  14-15  date.  As  discussed  previously, 
the  air  temperature  in  the  meadow  (Ml  site)  was  approximately  2°C 
lower  at  night  than  in  the  canopy  (M3  site)  during  two  July  dates. 

This  discrepancy  could  explain  some  of  the  deviation  between  simu- 
lated and  measured  results  during  the  night  period  in  Figure  12.  In 
addition,  as  discussed  previously,  it  is  believed  that  the  ERT's  of 
the  needle-branch  tips  may  bias  the  average  element  temperature  to  be 
high  which  would  also  explain  some  of  the  deviations  seen  in  Figure  12. 

Figures  12  and  13  show  that  the  simulated  temperatures  of  Layers 
1,  2,  and  3 closely  follow  air  temperature  during  the  day.  Unlike  the 
July  simulation.  Layer  1 temperature  does  not  deviate  significantly  from 
Layers  2 and  3 during  high  solar  irradiance.  This  is  due  to  the 
relatively  high  wind  velocity  and  thus  high  convectional  exchange  which 
occurs  during  the  October  simulation. 

To  recapitulate,  the  TCSM,  which  incorporates  the  geometric 
structure  of  a vegetation  canopy  and  predicts  the  thermal  response  of 
the  canopy  under  various  environmental  conditions,  was  developed. 

The  TCSM  is  designed  to  be  instantaneous  in  nature,  e.g.,  all  canopy 
elements  are  under  steady-state  conditions  and  no  heat  storage  may 
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Figure  12.  Simulated  versus  measured  lodgepole  pine  canopy 
temperatures  for  October  14-15,  1977.  Measured 
temperatures  are  the  mean  of  five  point  ERT's  of 
Layer  2 as  measured  by  the  Heat  Spy.  Air  temperature 
was  recorded  in  the  meadow  opening  (Ml  site) . 
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Figure  13.  Measured  global  solar  irradiance  (highest  value  of  Ml 
and  M2  sites)  and  air  temperature  (Ml  site)  for 
October  14-15,  1977. 
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occur.  Therefore,  the  model  is  independent  of  all  previous  environ- 
mental events.  In  many  applications  of  the  model  this  feature  would 
be  highly  desirable  (e.g.,  when  the  environmental  history  is  not 
known).  However,  to  approximate  the  time  dependent  phenomena  (non- 
steady-state  conditions),  one  can  use  a series  of  incremental  changes 
of  steady-state  energy  flows  (Gates,  1975).  In  branches  and  boles  a 
significant  amount  of  heat  storage  and  conduction  may  be  operating 
(Gates,  1975),  and  the  above  modifications  may  be  desirable  in  some 
applications. 

The  TCSM  was  applied  to  a lodgepole  pine  canopy.  The  algorithms 
incorporated  for  transpiration  and  convection  were  indeed  rather 
simplistic  in  their  assumptions.  For  example,  the  constant  R^ 
parameter  used  in  this  study  is  truly  variable  in  lodgepole  pine 
(Running,  197S).  One  of  the  greatest  barriers  in  applying  a model 
such  as  the  TCSM  to  a variety  of  vegetation  species  and  obtaining 
accurate  results  is  the  physiological  diversity  of  different  species, 
and  the  fact  that  the  physiological  response  of  many  species  are  not 
understood  sufficiently  to  be  predictable  (Running,  1978).  Depending 
on  the  researcher's  knowledge  of  the  vegetation  canopy  of  interest 
and  his  modeling  criteria,  more  appropriate  energy  transfer  algorithms 
can  be  incorporated  in  the  TCSM. 

In  the  future  the  model  can  be  made  more  comprehensive  in  nature 
by  including  algorithms  for:  thermal  radiant  reflections;  thermal 
spectral  radiance  of  the  sensor;  time  difference  equations  which 
account  for  energy  dynamic  and  heat  storage  of  the  soil  profile  and  the 
tree  boles;  wind  profiles  within  the  canopy;  air  temperature  profiles; 
and  canopy  water  relations.  In  addition,  it  is  anticipated  that  the 
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model  will  be  applied  to  a variety  of  vegetation  canopy  types  which 
do  not  closely  approximate  air  temperature  as  was  the  case  for  Lhe 
lodgepole  pine  canopy. 

Conclusions 

During  the  July  simulation  the  differentials  between  the  simu- 
lated versus  measured  horizontal  effective  radiant  temperatures  for 
a lodgepole  pine  canopy  were,  in  general,  less  than  2°C.  The  simu- 
lated average  layer  temperatures  for  all  three  layers  were  generally 
within  2°C  of  air  temperature,  except  for  Layer  1 during  periods  of 
high  solar  irradiance  in  July.  It  is  believed  this  discrepancy  is 
due  to  erroneous  values  of  the  internal  leaf  resistance  to  water 
vapor  diffusion  (R^).  The  simulated  resuits  suggest  that  for  needle- 
bearing  forest  canopies,  average  element  temperatures  deviate  signif- 
icantly above  air  temperature  only  during  periods  of  relatively 
high  solar  irradiance,  low  wind  velocities,  and  low  transpiration. 

The  effect  of  canopy  geometry  (element  inclination  distribution) 
on  thermal  radiation  transfers,  to  and  from  the  individual  layers 
within  the  canopy,  seems  to  be  minimal.  However,  the  canopy  geometry 
clearly  has  an  effect  on  the  contribution  of  thermal  radiation  from 
each  layer  to  the  sensor  above  the  canopy  as  a function  of  view  angle. 
It  is  believed  that  for  certain  canopy  element  inclination  distri- 
butions, canopy  LAI,  and  environmental  conditions,  the  sensor  incli- 
nation angle  will  affect  the  sensor  response  greatly;  and  this 
phenomenon  has  important  implications  on  the  optimum  view  angle  for 
making  inferences  about  the  target  of  interest. 
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The  surface  temperatures  of  canopy  elements  are  very  hetero- 
geneous, especially  during  direct  solar  irradiance  conditions.  The 
mean  of  several  point  measurements  (contact  thermistors  and  Heat  Spy 
on  branch  tips)  biases  the  measured  average  canopy  temperatures,  but 
a horizontal  ERT  which  includes  a cross  section  of  canopy  elements  is 
believed  to  be  the  least  biased  of  the  measurements  taken  in  this 
study. 

Both  air  temperature  and  wind  speed  variations,  as  a function 
of  location  within  the  canopy,  introduce  error  in  the  simulated 
results.  In  addition,  the  assumption  of  a constant  introduces 
error.  The  assumption  of  instantaneous  heat  exchange  may  not  be 
accurate  for  branches  and  boles. 

Due  to  the  large  heat  capacity  of  the  ground  and  net  thermal 
radiant  loss  to  the  clear  night  sky,  the  average  element  temperature 
increases  as  one  proceeds  from  Layer  1 to  Layer  3 as  demonstrated  by 
both  simulated  and  measured  results. 

Model  simulations  showed  that  the  total  global  irradiance  ab- 
sorbed by  the  lodgepole  pine  canopy  system  is  relatively  constant  with 
solar  zenith  angle.  However,  the  proportion  of  total  global  irradi- 
ance absorbed  by  individual  canopy  layers  varies  as  a function  of 
solar  zenith  angle. 

The  TCSM  provides  a modeling  framework  which  may  be  useful  to  a 
variety  of  research  interests.  Specific  energy  transfer  algorithms, 
which  are  best  suited  to  the  researcher's  modeling  criteria,  can  be 
incorporated  in  the  TCSM.  The  TCSM  is  unique  in  that  the  framework 
incorporates  the  geometric  structure  in  radiation  transfer  algorithms. 
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4.0  CONCLUSIONS  AND  RECOMMENDATIONS 


■t 


The  thermal  canopy  signature  model  (TCSM) , Chapter  3,  together 
with  the  absorption  calculations  performed  by  A15S0RPT,  Chapter  2, 
provide  a modeling  framework  for  simulating  energy  transfers  within 
vegetation  canopies  of  specific  geometric  structure.  The  thermal 
modeling  approach  is  unique  in  that  it  incorporates  detailed  geometric 
canopy  structure  to  define  the  radiant  transfers  occurring  within  the 
canopy  system.  The  model  predicts  the  radiometric  temperatures  as  a 
function  of  look  angle.  In  addition,  a radiometric  temperature  height 
profile  of  the  canopy  is  calculated.  The  model  was  successfully 
applied  to  a lodgepole  pine  (Pinus  contorts)  canopy  which  indicated 
that,  except  for  specific  environmental  and  physiological  conditions, 
mean  element  temperatures  closely  approximate  air  temperature.  The 
model  suggests,  however,  that  canopy  structure  can  have  significant 
effects  on  the  response  of  a thermal  sensor  above  the  canopy.  The 
angular  thermal  exitance  prediction  capability  of  the  TCSM  recommends 
it  as  a useful  tool  for  defining  optimum  sensor  view  angle  and  environ- 
mental conditions  for  target/background  discrimination  studies. 

The  research  described  in  this  report  also  addressed  the  devel- 
opment of  methods  for  determining  forest  canopy  geometry,  particularly 
needle  and  branch  angle  frequency  distributions.  It  was  found  that 
these  distributions  strongly  influence  the  absorption  of  solar  radiation 
within  the  canopy.  Solar  radiant  absorption,  in  turn,  affects  the 
overall  thermal  equilibrium  within  the  canopy  layers.  A description 
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of  the  geometry  techniques  lias  been  submitted  and  accepted  by  the 
appropriate  journals  (Appendix  A). 

There  are  many  broad  areas  of  potential  application  of  the 
modeling  approaches  described  depending  on  a user's  orientation.  Some 
of  these  applications  might  focus  on  the  vegetation  canopy  itself  as 
the  primary  scene  of  interest,  as  in  agricultural  or  forest  water 
relation  studies.  Others  might  focus  on  the  surface  beneath  the 
canopy,  as  in  geological  or  snow  cover  estimation  problems.  A detailed 
exploration  and  description  of  such  applications  is  beyond  the  scope 
of  this  report.  However,  with  regard  to  the  modeling  approach  itself, 
the  authors  make  the  following  recommendations. 

First,  and  most  importantly,  how  well  does  the  TCSM  perform  when 
applied  to  a wide  variety  of  vegetation  types  and  terrain  conditions? 

The  model  needs  to  be  validated  for  other  situations  than  the  lodgepole 
canopy  simulated  here.  Undoubtedly , specific  energy  transfer  algorithms 
presently  incorporated  in  the  model  will  need  modification.  The 
application  of  the  model  to  closed  versus  open  canopies  will  also 
need  to  be  systematically  addressed. 

Secondly,  what  is  the  tradeoff  in  model  precision  or  accuracy 
versus  the  availability  of  model  parameters ? A related  question  is 
the  relative  advantage,  under  different  conditions,  of  employing  a 
detailed  model  such  as  the  TCSM  which  meticulously  accounts  for  all 
energy  transfers  and  incorporates  complete  canopy  geometry  to  that 
of  employing  more  simple  models? 

In  summary,  the  terrain  feature  models  described  in  this  report 
provide  a basic  framework  which  can  be  adapted  or  utilized  to  study 
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a large  number  of  research  interests.  For  many  design  problems 
relative  to  predicting  background  electromagnetic  behavior  of  natural 
features,  the  present  models  should  prove  useful.  However,  further 
validation  and  appropriate  modification  of  model  processes  is 
recommended  before  the  model  is  applied  generally. 


A Monte  Carlo  Calculation  of  the  Effects  of 
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ABSTRACT 

The  ability  of  vegetation  canopies  to  absorb  photosynthetically 
active  radiation  (PhAR)  is  known  to  be  a function  of  the  canopy  geometry. 
A Monte  Carlo  model  was  used  to  estimate  relative  PhAR  absorption  in 
theoretical  vegetation  canopies  of  different  structure.  The  relatively 
simple  single  component,  multilayer  simulations  adequately  describe  the 
empirically  established  trends  reported  by  a wide  variety  of  investi- 
gators. Absorption  trends  for  erectophile  (mostly  vertical  leaves)  and 
planophile  (mostly  horizontal  leaves)  canopies  are  indicated  with 
respect  to  leaf  area  index  (LAI)  and  solar  zenith  angle.  Generally, 
our  model  results  show  that  erectophile  canopies  are  more  efficient  at 
absorbing  PhAR  under  medium  to  high  LAI  and  all  ranges  of  zenith  angle. 
Planophile  canopies  show  increased  absorption  at  lower  LAI's. 
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Abstract 

Optical  diffraction  analysis  of  i_n  situ  ground  photographs  has 
previously  been  utilized  to  estimate  foliage  angle  distributions  in 
grassland  canopies.  These  canopies  are  typically  characterized  by  a 
single  component,  leaves,  and  the  foliage  is  highly  linear  in  nature. 

In  this  paper,  the  diffraction  technique  is  extended  to  a r ilti-component 
forest  canopy  containing  needles  and  branches.  Additional  convolution 
and  coordinate  transformations  are  derived  to  estimate  the  branch  and 
needle  angle  frequency  distributions  for  top,  middle,  and  base  sections 
of  two  lodgepole  pine  (Pinus  contorta)  trees.  The  resulting  distribu- 
tions show  that  the  branch  inclination  angles  tend  to  increase  as  one 
proceeds  to  the  tree  tops.  The  needle  inclination  angle  distribution 
was  relatively  constant  for  all  layers,  and  it  is  believed  that  this 
distribution  is  characteristic  of  a large  class  of  needle  bearing 
species . 
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Abstract 

A non-destructive,  rapid  technique  utilizing  horizontal  ;n  sun  ground  photographs  for  estimating 
foliage  angie  distributions  :s  discussed.  Optical  diffraction  patterns  generated  from  orthogonal 
photographs  are  analysed  for  ar.guur  mas  oy  wedge  sampling.  Probability  distributions  for  planar 
projections  of  foliage  orientations  are  derived  from  these  measurements  and  mathematically  convo- 
luted to  determine  the  actual  three-space  propability  distribution  function  for  foliage  angles.  The 
method  is  particularly  appropriate  for  dense  canopies  which  are  difficult  to  measure  by  other  tech- 
niques. The  diffraction  technicue  :s  evaluated  for  abstract  canopies  and  for  a canopy  of  Western 
wheat  grass  ( Agropvron  smuhti).  It  also  yielcs  physically  consistent  nterpretations  for  the  phenolo- 
gies! development  of  domestic  Satanta  wheat  ( Tnncitm  aertivnm). 


Introduction 

Foliage  angle  distribution  functions  for  grassland  canopies  have  traditionally  been 
estimated  by  point  quadrat  techniques  such  as  discussed  by  Wilson  ( I960.  1963)  and 
Philip  11965).  Other  recent  techniques  involving  photographic  or  photocell  measure- 
ments of  foliage  gap  frequency  which  can  be  related  to  the  foliage  angle  distribution 
include  the  methods  discussed  by  Norman  and  Tanner  { I C)69')  and  Bonhomme  and 
Chartier  (1972).  The  present  authors  described  another  photographic  procedure  for 
grassland  canopies  in  an  earlier  paper  i Smith  ei  ai.  1977)  whereby  orT-angie  photographs 
are  used  to  record  gap  frequency  and  a Fredholm  integral  equation  is  solved  which 
relates  foliage  angles  to  gap  frequency. 

The  point  quadrat  technique  offers  practical  difficulties  in  the  amount  of  held  time 
required  to  obtain  the  measurements.  The  photographic  and  photocell  methods  are 
an  improvement  in  this  regard,  but  they  are  difficult  to  apply  in  canopies  with  dense 
foliage  cover.  The  Fredholm  technique,  for  example,  cannot  be  applied  to  dense 
canopies  in  which  canopy  closure  at  most  view  angles  is  complete. 

In  this  paper  we  present  an  alternative  approach  applicable  to  dense  canopies 
which  utilizes  optical  diffraction  pattern  analysis  of  held  photographs.  Planar  distri- 
butions of  foliage  angles  are  determined  from  orthogonal  ground  photographs  obtained 
in  the  vertical  plane  of  the  plant  canopy.  These  orthogonal  distributions  are  then 
mathematically  convoluted  to  estimate  the  actual  plant  canopy  foliage  angie  distri- 
bution. This  diffraction  technique  is  also  a simple  and  rapid  in  situ  measurement 
method. 
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Abstract 

There  ts  an  increasing  interest  in  theoretical  models  which  describe  the  interaction  of  solar  radiation 
with  vegetation  canopies.  Common  to  these  models  is  a need  to  describe  mathematically  the  geo- 
metric structure  of  the  plant  canopy  The  amount  of  radiation  reflected  or  absorbed  by  the  canopy 
is  primarily  determined  by  the  distribution  of  gaps  in  the  foliage  with  respect  to  the  radiation  source. 
A measure  of  canopy  geometry  related  to  gap  frequency  at  various  view  angles  is  the  distribution  of 
leaf  angles.  Two  methods  for  measuring  the  distribution  of  leaf  angles  are  discussed.  The  first 
method  is  to  project  orthogonally  and  photograph  individual  plants  and  relate  the  measured  leaf 
angles  in  the  projections  to  the  canopy  distribution  of  angles.  The  second  method  is  a rapid  in  situ 
method  based  on  ground  level  multiple  view  angle  photography.  A Fredholm  integral  equation 
relating  foliage  angles  to  the  proportion  of  gap  in  the  canopy  as  a function  of  view  angle  is  then 
solved.  Comparisons  of  the  results  using  the  two  methods  are  made  for  a canopy  of  Western  wheat 
grass  ( Aaropvron  smithii). 


Introduction 

Mathematically,  a homogeneous  vegetation  canopy  may  be  described  given  the 
following  information: 

(1)  inclination  angle  distribution  of  the  foliage  elements. 

(2)  Azimuthal  angle  distribution  of  the  foliage  elements. 

(3)  Leaf  area  index  (ratio  of  the  one-sided  leaf  area  to  a unit  area  of  underlying 
soil  surface). 

(4)  A relation  describing  the  three-dimensional  dispersion  of  the  leaf  area  within 
the  canopy. 

The  situation  becomes  more  complex  if  the  canopy  is  heterogeneous  in  either 
composition  or  structure.  Heterogeneity  usually  implies  that  the  canopy  must  be 
stratified  into  layers  and  the  above  information  determined  for  each  layer.  Stratifica- 
tion may  be  determined  either  from  a height  distribution  (Oliver  and  Smith  197?)  or 
from  the  apparent  morphology  characteristics  of  the  vegetation  under  study. 

The  foliage  inclination  angle  distributions  for  various  typical  classes  of  stand 
geometry  arc  shown  in  Fig.  I (dc  Wit  1965).  Horizontal  leaves  are  most  frequent  in 
planophilc  canopies,  anti  vertical  leaves  in  crcctophile  canopies.  The  leaves  in 
plagiophilc  canopies  arc  most  frequent  at  oblique  inclinations  of  greater  than  45  , 
* those  in  extremophile  canopies  at  oblique  inclinations  of  less  than  45  . Interpretations 

of  the  cumulative  frequency  distribution  function  may  be  made  as  above  and  their 
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Abstract 

A battery-powered  scene  recording  radiometer  system  has  been 
developed  for  relating  spectral  variability  and  target  composition. 
A remote  controlled  filter  wheel  radiometer  is  interfaced  with  a 
llasselblad  500  EL  camera,  so  that  the  silicon  detector  is  directed 
toward  the  camera's  viewing  glass  Signal  detection  at  discrete 
wavelength  bands  is  achieved  by  successively  rotating  interchange- 
able interference  filters  that  interdict  the  view  of  the  detector. 
Filter  positioning  is  controlled  by  coding  holes  drilled  in  the  fil- 
ter wheel  disk  which  are  interpreted  by  a bank  of  opposing  light 
emitting  diodes  and  phototransistors.  Upon  obtaining  a photo- 
graphic record  of  the  scene,  the  camera  is  automatically  advanced. 

Introduction 

A recurring  problem  in  applying  remote  sensing  technology  to 
natural  resources  is  the  difficulty  in  correctly  classifying  a reso- 
lution element  which  contains  a mixture  of  materials.  This  prob- 
lem is  particularly  acute  in  natural  vegetation  communities  such 
as  those  encountered  in  the  West,  where  a great  deal  of  hetero- 
geneity occurs  within  the  resolution  clement  A machine-assisted 
classification  rule,  however,  is  forced  into  one  of  two  choices: 
(1)  identifying  the  resolution  element  as  a single  material  when, 
in  fact,  it  may  contain  only  a small  percentage  of  the  material, 
or  (2)  leaving  the  resolution  clement  unclassified.  This  problem 
of  mixtures  is  closely  related  mathematically  to  the  problem  of 
signature  extraction,  which  describes  the  difficulty  of  determin- 
ing a typical  spectral  response  for  materials  when  the  underlying 
data  distributions  are  heterogeneous. 

There  are  two  broad  approaches  to  the  mixtures  problem. 
These  include  a least  squares  approach,  and  a parameter  esti- 
mation technique  using  maximum  likelihood  procedures.  Pace 
and  Detchmendy1  as  well  as  Hallum2  discuss  the  former  ap- 
proach. The  maximum  likelihood  procedure  is  described  by 
Horwitz  et  al3  and  Guseman".  Both  approaches  have  as  their 
fundamental  assumption  the  hypothesis  that  there  is  a linear 
relationship  between  measured  spectral  response  from  a reso- 
lution element  and  the  proportions  of  materials  contained  with- 
in the  resolution  element  This  assumption  has  never  been  clearly 
evaluated  through  empirical  data  The  primary  reason  for  this 
omission  appears  to  be  the  fact  that  most  spectrometers  measure 
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the  total  response  from  a resolution  element  without  the  capa- 
bility of  providing  a registered  record  of  the  composition  of  the 
resolution  element.  One  exception  is  the  system  reported  by 
Dana5  as  used  in  determining  aircraft  reflectance  measurements 
in  which  a bore-sighted  camera  is  used  in  conjunction  with  a 
radiometer. 

The  instrument  described  in  this  paper  was  developed  to  pro- 
vide a method  for  obtaining  a direct  and  simultaneous  measure- 
ment of  both  scene  composition  and  spectral  response,  using  a 
common  field  of  view.  The  instrument  has  been  used  by  the  au- 
thors for  investigating  mixture  effects  in  a lodgcpole  pine  ( I'inus 
corttorta)  community  in  the  central  Colorado  mountains.  The  age 
of  the  stand  studied  was  between  30  and  50  years  with  average 
height  of  the  stand  approximately  20  ft.  Canopy  density  was  van- 
able  ranging  from  approximately  8057  crown  closure  to  a com- 
pletely open,  grass-covered  clearing.  The  majority  of  the  under- 
story within  denser  regions  of  the  lodgepole  stand  was  similar 
to  that  of  the  meadow  opening  In  this  application  an  aerial  tram- 
way system  was  required  to  suspend  the  instrument. 

Physical  Features 

The  scene  recording  radiometer  is  essentially  a remote  controlled 
filter  wheel  radiometer  interfaced  to  a Hasselblad  EL  500  camera. 
A schematic  of  the  instrument  is  shown  in  Figure  1.  A remote 
readout  and  control  station  are  utilized  in  conjunction  with  the 
instrument.  The  functions  of  the  system  arc  five  fold  ( 1 ) signal 
detection,  (2)  partitioning  of  the  signal  into  discrete  wavelength 
bands,  (3)  control  of  filter  wheel  position,  (4)  photographing  the 
scene,  and  (5)  remote  system  control  and  signal  readout. 

The  radiometer  system  uses  a silicon  detector  directed  toward 
the  viewing  glass  of  the  camera.  The  optics  of  the  camera  define 
the  instrument's  field  of  view.  This  allows  adjustment  of  the  ra- 
diometer field  of  view  by  varying  the  camera's  lens  system.  Sharp 
definition  of  the  field  of  view  of  the  instrument  is  achieved  by 
overlaying  a circular  mask  on  the  viewing  glass.  The  diameter  of 
the  mask  was  constructed  to  be  slightly  smaller  than  the  physi- 
cally constrained  field  of  view  of  the  detector  when  interdicted 
by  an  interference  filter.  This  constrained  field  of  view  was  mea- 
sured by  mounting  the  instrument  on  a goniometer  and  rotating 
it  about  the  optical  axis  of  the  camera,  while  exposing  it  to  a 
narrow  beam  of  light.  An  mtcrference  filter  centered  at  6328  A 
was  used  to  accommodate  a helium  neon  gas  laser  light  source. 
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Scene  Radiation  Dynamics 


ABSTRACT 


These  reports  are  a two-volume  Final  Report  Series  for  Project 
DACW  39-77-C-0073  issued  by  the  Environmental  Laboratory  of  the  U.S. 

Army  Engineer  Waterways  Experiment  Station.  The  period  covered  is  from 
11  August,  1977  to  30  September,  1978.  Overall  objectives  of  this  one 
year  study  were  to  develop  a comprehensive  optical  and  thermal  signature 
data  base  and  to  evaluate  or  develop  optical  and  thermal  canopy  radia- 
tion models.  A variety  of  vegetation  terrain  features  were  studied 
including  coniferous  trees  (Pinus  c ontorta) , deciduous  trees  (Elaeagnus 
augustifolia) , shrubs  (Potentilla)  and  grasses  (Festuca) . In  order  to 
synthesize  the  scene  radiation  dynamics  with  the  use  of  models,  accom- 
panying geometric  and  meteorological  parameters  were  also  obtained. 

Volume  I presents  the  optical  and  thermal  modeling  descriptions  and 
includes  terrain  data  modules  for  coniferous,  deciduous,  and  grass  ter- 
rain features.  The  optical  SRVC  (Solar  Radiation  Vegetation  Canopy 
Model)  developed  under  previous  U.S.  Army  Research  Office  sponsorship 
is  evaluated  against  these  terrain  modules.  A thermal  leaf  model  and 
an  initial  thermal  canopy  signature  model  are  described  and  compared 
against  field  measurements.  Both  optical  and  thermal  signature  models 
are  infinite  plane  terrain  approximations  to  a three-layer  stratified 
canopy.  Source  and  view  angle  dependencies  of  the  exitance  are  pre- 
dicted. In  addition,  the  thermal  model  predicts  the  temperature  distri- 
bution of  the  vegetation  layers. 

Volume  II  contains  the  optical  reflectance  data  listings  and  in- 
cludes descriptive  information  for  the  experimental  sites.  The  data 
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types  and  data  reduction  methods  are  enumerated.  Finally,  an  analysis 
of  optical  data  dispersions  is  given  including  seasonal  and  diurnal 
variability  and  two-spectral  space  scatter  plots.  The  applicability  of 
a Tasseled  Cap  type  of  transformation  is  evident. 
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Evaluation  of  Illumination  and  Terrain  Geometry  Effects  on 
Spectral  Response  in  Mountain  Terrain 

ABSTRACT 

An  extensive  analysis  of  terrain  geometric  effects  on  the  optical 
scattering  properties  of  natural  resource  scene  in  mountainous  terrain 
has  been  performed.  Spectral  reflectance  measurements  were  obtained 
for  lodgepole  pine,  Pinus  contorta,  ponderosa  pine.  Firms  ponderosa, 
Russian  olive,  Eleaegnus  angusti folia,  grass  species,  Agvopyron  sp. , 
and  Eronrus  so. , and  snow.  Sensor  platforms  included  ground-based 
measurements  using  aerial  tramways,  aircraft  radiometric  observations, 
and  satellite  (Landsat)  measurements.  A wide  range  of  effective  view 
and  source  illumination  angles  were  recorded  for  the  various  target/ 
sensor  combinations. 

Regression  analyses  and  photometric  plots  were  made  from  the  data 
in  order  to  test  the  Lambertian  assumption  for  the  various  material 
types.  In  addition  a process-oriented  radiative  transfer  model  was 
applied  to  the  data.  This  model  was  also  used  to  evaluate  initial 
effects  of  background  topographic  variations. 

Results  of  this  study  indicate  that,  particularly  in  the 
chlorophyll  absorption  band  all  materials  exhibit  non-Lambertian 
behavior  for  effective  zenith  sensor  or  source  angles  greater  than 
60  degrees,  but  that  for  effective  angles  less  than  AO  degrees,  the 
Lambertian  assumption  may  be  valid.  For  stable  atmospheric  conditions 
and  constant  phase  angle  the  Minnaert  relationship  may  be  applied  to 
quantify  scene  radiance  properties.  The  canopy  reflectance  mode]  was 
found  to  follow  the  general  trends  of  the  field  measurements  but  over- 
estimates infrared  response.  In  order  to  adequately  model  topographic 
influences  or  spectral  response,  canopy  density  variations  must  be 
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APPENDIX  B:  Supporting  Material  for  Monte  Carlo  Calculations  of 
Absorbed  Solar  Radiation  in  Vegetation  Canopies 


Diffuse/Direct  Ratio  Sensitivity  Analysis 
Element  Transmission  Sensitivity  Analysis 
Vertical  Reflection  Validation 
SRVC  Simulated  Absorption  Coefficients 


Program  Listing  for  ABSORPT 


APPENDIX  B 


The  following  topics  present  in  f ill  detail  the  results  of  the 
data  analysts  and  computer  programs  mentioned  in  the  main  text  of 
Chapter  III. 

Dif fuse/Direct  Pa tjx  Sen s_L tivitv  Analysis 

Figures  1 and  2 present  the  proportion  of  absorbed  spectral  solar 
irradiance  for  a wavelength  of  high  canopy  reflectance  (p  = .50, 

T = .05  for  all  canopy  elements)  with  a solar  zenith  angle  of  0°  and 
12°,  respectively.  Figures  3 and  4 present  the  proportion  of  absorbed 
spectral  solar  irradiance  for  a low  canopy  reflectance  (p  = .05, 
i = .001  for  all  canopy  elements)  with  a solar  zenith  angle  of  0°  and 
12°,  respectively. 

Element  Transmission  Sensitivity  Analysis 

Figures  5 and  6 show  the  simulated  proportion  of  absorbed  spectral 
solar  irradiance  for  a wavelength  of  high  canopy  reflectance  (p  = .50 
for  all  canopy  elements)  as  a function  of  canopy  element  transmission 
coefficient  for  a solar  zenith  angle  of  0°  and  12°,  respectively. 

Vertical  Reflection  Vnl i da t i on 

The  measured  versus  simulated  vertical  spectral  canopy  reflectance 
of  the  modeling  trees  for  three  time  periods  are  presented  in  Figures 
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S RVC  f irm  1 a t yd  Abs  orjU  ion  Coefficients 

The  Simulated  . absorption  coefficients  for  various  mean 
canopy  element  reflectances  and  for  a solar  angle  of  4 7 are  shown 
in  Figure  10. 

P r o i r a n ABSORPT 

The  Fortran  program  ABSORPT  calculates  the  spectral  and  total 
absorption  within  the  canopy  system  using  the  interpolation  and 
integration  techniques  described  in  Chapter  LIT. 


Absorption  Coefficient 


Hif fust /Direct  Ratio 


Figure  1.  The  simulated  proportion  of  absorbed  spectral  solar 

irradiance  for  a wavelength  of  high  canopy  reflectance 
(p  = .50,  x = .05  for  all  canopy  elements)  with  a sola: 
zenith  angle  (/.)  of  0°.  The  proportion  absorbed  (a.) 
by  each  canopy  layer  is  shown  as  well  as  the  total  1 
canopy  absorption. 


Diffuse/Direct  K.ilio 


Figure  3.  The  simulated  proportion  of  absorbed  spectral  solar 

irradiance  for  a wavelength  of  low  canopy  reflectance 
(p  = .05,  t = .001  for  all  canopy  elements)  with  a 
solar  zenith  angle  of  0°.  The  proportion  absorbed  (a.) 
by  each  canopy  layer  is  shown  as  well  as  the  total 
canopy  absorption. 


Absorption  Coefficient 


Figure  4.  The  simulated  proportion  o£  absorbed  spectral  solar 

irradiance  Cor  a wavelength  of  low  canopy  reflectance 
(p  = .05,  t = .001  for  all  canopy  elements)  with  a 
solar  zenith  angle  of  72°.  The  proportion  absorbed  (c<.) 
by  each  canopy  layer  is  shown  as  well  as  the  total  1 
canopy  absorption. 


Absorption  Coefficient 


Transmission  of  Elements 


Figure  5.  The  simulated  proportion  of  absorbed  spectral  solar 

irradiance  for  a wavelength  of  high  canopy  reflectance 
(p  = .50  for  all  canopy  elements)  for  four  canopy 
element  transmission  coefficients  and  for  a solar  zenith 
angle  of  0°.  The  proportion  absorbed  (a.)  for  each  canopy 
layer  is  shown  as  well  as  the  total  canopy  absorption. 


Figure  6. 


Transr.ir.sion  of  Elements 


The  simulated  proportion  of  absorbed  spectral  solar 
irradiance  for  a wavelength  of  high  canopy  reflectance 
(p  = .50  for  all  canopy  elements)  for  four  canopy  element 
transmission  coefficients  and  for  a solar  zenith  angle 
of  12°.  The  proportion  absorbed  (a.)  for  each  canopy 
layer  is  shown  as  well  as  the  total  canopy  absorption. 


Kef lectnnce 


Wavelength  (v»w) 


Figure  8.  Measured  versus  simulated  vertical  spectral  canopy 
reflectance  for  July  1,  1977,  1000  Standard  Time. 


Mean  Canopy  Element  Reflectance  (MCR) 


Figure  10.  SRVC  simulated  a.  . absorption  coefficients  for  various 

A , 1 , Z 

ir, can  canopy  element  reflectances  and  for  a solar  zenith 

angle  of  47°.  The  a , for  each  canopy  layer  is  presented 

A , l > z 

as  well  as  the  total  proportion  of  the  spectral  irradiance 
absorbed . 
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APPENDIX  C 


The  following  topics  present  in  full  detail 
and  computer  programs  mentioned  in  the  main  text 


the  data,  analyses 
of  Chapter  IV. 
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Data  Reduction  and  Initial  Ar ■ .lysis 

Figures  1 and  2 present  the  measurements  of  air  temperature 
(Ml  Site),  the  contact  thermistor  and  the  Heat  Spy  on  two  branch  tips 
for  October  14-15,  1977.  Figure  3 shows  the  simulated  TCSM  results 
for  October  14-15  versus  the  three  contact  thermister  measurements. 

Dav  and  Night  Input  and  Output 

The  TCSM  input  and  output  for  0930  and  0330  Standard  Times  for 
July  15-16,  1977,  are  presented  in  Figures  4 and  5.  The  input  data 
for  the  day  and  night  environmental  conditions  are  the  fixed 
parameters  used  for  the  sensitivity  analysis. 

Sensitivity  Analysis 

Figures  6,  7,  8,  and  9 show  the  emissivity  and  internal  resistance 
to  water  vapor  diffusion  sensitivity  results  for  the  day  and  night 
environmental  conditions. 

Air  Tempera ture  Variations 

Figures  10  and  11  show  the  air  temperature  measurements  for  Ml, 

M2,  and  M3  Sites  for  July  14  and  July  15-16,  1977.  Figure  12  presents 
the  simulated  versus  measured  horizontal  ERT's  for  July  15-16.  The 
air  temperature  was  recorded  in  the  meadow  opening  (Ml  Site)  and  is 
presented  in  Figure  13. 
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Program  TC'SM 

The  Fortran  program  TCSM  is  the  thermal  canopy  signature  model, 
as  presented  in  Chapter  IV.  The  required  inputs  are  described  in 
Fubroutine  Inputda.  The  program  ZSVSTM  of  the  International  Mathe- 
matical and  Statistical  Library  (1977)  must  be  attached. 


6 11  16  21  2 7 12 

Standard  Time 


Figure  1.  Measured  contact  thermister  and  Heat  Spy  data  for  a branch 
tip  in  Layer  2,  and  air  temperature  at  the  Ml  Site  for 
October  14-15,  1977. 
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Figure  3.  Simulated  mean  canopy  layer  temperature  versus  the  data 
from  three  contact  thermistors  in  branch  tips  of  Layers 
1,  2,  and  3 for  October  14-15,  1977.  Air  temperature  was 
recorded  at  the  Ml  Site. 
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* Time  = 0930 

* Eniissivity  (Layers  1-4)  = 1.0,  1.0,  1.0,  1.0 

2 

* Solar  irradiance  = 855  w/m 

* Air  temperature  = 14.6°C 
Ground  temperature  = 11.7°C 

* Wind  velocity  = 10.0  cm/s 

* Relative  humidity  = .20 

* Leaf  resistance  to  water  vapor  diffusion  = 0.66  nin/cm 

Average  element  temperatures  (Layers  1-3)  = 20.4,  16.6,  16.3°C 

2 

Ground  thermal  exitance  = 372  w/m 

2 

Sky  thermal  exitance  = 302  w/m 

* Mean  absorbed  solar  flux  density  for  elements  Layers  1-3  = 

144,  49,  46  w/m~ 

Emitted  thermal  flux  density  for  the  average  elements  (1-3)  = 
420,  399,  397  w/m2 

Absorbed  thermal  flux  density  for  the  average  elements  (1-3)  = 
375,  387,  384  w/m2 

Energy  gain  by  convection  for  the  average  elements  (1-3)  = 

-90,  -31,  -26  w/m1- 

Energy  loss  by  transpiration  for  the  average  elements  (1-3)  = 
9,  7,  7 w/m2 


The  thermal  exitance  and  ERT 

viewing  angles  are: 

Inclination  (degrees) 

above  the  canopy  for  the 
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The  thermal  exitance  and  ERT 

for  the  horizontal  view 

are: 

Layer 

Exitance  (w/m 

ERT  (°C) 

1 

420 

20.4 

2 

399 

16.6 

3 

397 

16.3 

Figure  4.  Environmental  and  simulated  conditions  for  0930  (day  environ- 
mental conditions)  July  15,  1977.  Input  parameters  are 
denoted  by  asterisks. 
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Time  = 0330 

Emissivity  (Layers  1-4)  - 1.0,  1.0,  1.0,  1.0 

2 

Solar  irradiance  «*  0.0  w/m 
Air  temperature  = 0.4°C 
Ground  temperature  = 5.0°C 
Kind  velocity  = 10  cm/s 
Relative  humidity  = 0.S5 

Leaf  resistance  to  water  vapor  diffusion  = 0.66  min/cm 

Average  element  temperatures  (Layers  1-3)  = -1.0,  0.0,  0.5°C 

2 

Ground  thermal  exitance  = 399  w/m 

2 

Sky  thermal  exitance  = 234  w/m 

2 

Mean  absorbed  solar  flux  for  Layers  1-3  = 0.0,  0.0,  0.0  w/m 
Emitted  thermal  flux  density  for  the  average  elements  (1-3)  = 
310,  315,  317  w/m2 

Absorbed  thermal  flux  density  for  the  average  elements  (1-3)  = 
288,  309,  320  w/m2 

Energy  gain  by  convection  for  the  average  elements  (1-3)  = 

22,  7,  -2  w/m2 

Energy  loss  by  transpiration  for  the  average  elements  (1-3)  = 
0.2,  0.4,  0.5  w/m2 


The  thermal  exitance  and  ERT  above  the  canopy  for  the  various 
viewing  angles  are: 


ion  (de 

grees) 

Exitance  (w 

hr) 

ERT 

(°C) 
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-0 
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.1 

55 
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-0 

.1 

65 
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-0 

.1 

75 

315 

-0 

.1 

35 

315 

-0 

.1 

:xitance 

and  ERT 

for  the  horizontal 

view  are: 

Layer 

Exitance  (w 

7m2) 

ERT 

(°C) 

1 

310 

-1 

.0 

2 

315 

0 

.0 

3 

317 

0 

.5 

igure  5.  Environmental  and  simulated  conditions  for  0330  July  16, 
1977.  Input  parameters  are  denoted  by  asterisks. 
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Figure  6.  Sensitivity  analysis  of  average  element  emissivity  versus 
average  clement  temperature  of  the  three  canopy  layers 
for  the  day  environmental  conditions. 
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Figure  7.  Sensitivity  analysis  of  average  element  cmissivity  versus 
average  element  temperature  of  the  three  canopy  layers 
for  the  night  environmental  conditions. 


onperaiure 


o .2  .4  .6  .8  1.0  1 . 2 


Resistance  to  Water  Vapor  Diffusion  (run, 'em) 


Figure  8.  Sensitivity  nnalysis  of  average  element  resistance  to 

water  vapor  diffusion  versus  average  element  temperature 
of  tlie  three  canopy  layers  for  the  day  environmental 
cond it  ions . 
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Figure  9.  Sensitivity  analysis  o£  average  element  resistance  to 

water  vapor  diffusion  versus  average  element  temperature 
of  the  three  canopy  layers  for  the  night  environmental 
conditions . 


Figure  10.  Air  temperature  measurements  at  the  Ml,  M2,  and  M3  Sites 
for  July  14,  1977. 
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Figure  11.  Air  temperature  measurements  at  the  Ml,  M2,  and  M3  Sites 
for  July  16,  1977. 
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Figure  12.  Simulated  versus  measured  lodgepole  pine  canopy  horizontal 
effective  radiant  temperatures  (KRT)  for  July  15-16,  1977. 
Measured  F.RT's  are  the  mean  of  4 horizontal  ERT's  of  the 
middle  layer  as  measured  by  the  Heat  Spy.  Air  temperature 
was  recorded  at  the  Ml  Site. 
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GO  TO  130 

120  CONVEC*  I 0. 9S  *wv*  *0 . 97  ) • I . 00  1 ) • I 69  7.  7M 


4iiii  IS  tSkST  QUALITY  KhACIlUk*U* 

iHOU  OJrX  FMklSHjCL  TO  UDC  . 
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CnNN/f  C • CONMf ( -1  . 1 
c:»Nvr  c »c  cnvc:  • < x c*<  i -ah 
130  CONTINUE 

65  r 

C % . . F UNC T I ON  TMMLEX  CALCULATES  T m F t**TTTCf>  Tmto-al  FNfBGY  fBO"  THE  AVftAf.t  LfAf 
r IN  L 'Tf  3 K. 

c 

H^LFXi  $TEe#SMISSV<“»MM*l*N»**A 

70  C 

C... FUNCTION  T»AN*  r ALCiiL  * TE  S Thf  -NrcoY  LOSS  F3Hm  Tug  A Vf  » AGE  LEAF  IN  lAVEB  it 
c 5 Y T^ANSPIft  ATTTN. 

C 

l f • .»0 . 5f  5 • V ( *1  » KT7.3 

75  ^ )L  ■ C 5 . ? 3N’*  F * 0 (O.O^4' 715  • v (<  M ) •1.0F-6 

S 0 A ■ O‘4*(5.',X4.7*F»r(C.05^715«AT)  1 M . 0 1 
3 a • I0.0*M.’7M  1./wV**T.sn/bC. 

TwANc»Lf*(5«t»r0A»/(OL*®AI 

T»  ANS-6W.75  • T»i«^ 

80  C 

C...SUN  AIL  THF  6N=BGY  L0S5f5  AND  r.  A I N * OP  Th:  A V ? B A G E LEAF  IN  L A Y £ 3 K « 

c 

TOTAL-  THEa"-THMtcy*mNV£:-T3ANS*AaSnL  ( « ) 

3-  TUBN 

05  t NO 


,1* 


. r>v' 


C^vv; 


c i 

,4»V\  V ^ • 

- jf  v 

vl  -•  ^ 
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surpputine  g r viiri 
c 
c 

C . . . SUAQOUI I Hf  OVA  *l  G CAt^i'LAT-S  TH  c Coc  ( L F > Oc  V I A T I ON 

C as  Gl  E OF  AIL  L c A r isri  IS4T  m • p *-  r POIFhtmqn<  r[  P“UT  A T I P^S  . Inf  THE  QPY 

C IS  RASED  ON  Thc  FXISTFNCF  nc  pi  nc  FlE-PNTS  AS  USED  IN  THE  SPVC  MODEL. 


COMMON  /a/  UV  , PH.  ® L . T f R 1 

COM  - ON  / R / PSAP(1.9»,  OHirn,Q),PraP?n.O),  PMlT?r»,Q).STEF 
COM-ON  /c /:  OS  TA  ( Q,  9,  } « ) , - M ! s c V I A | , ARS  °P  R ( ?)  . ES<y,  E Gc  0 » SFCTap(Q) 
common /o /c  p*iT  ( ?•  *.  0 ( 1 . 5 • o ) , Sl>MTf3t9)»  xtLVi  GT,  ITIf*F 

CrMMON/F<AT»THFTA(<5).PHt(lfl>«*LP('>).  VLF(9),  7 L ? « «*  I * * S ( 9 , 1 R ) 

C n M.  m pn  / F / T $ ( ? * 1 R ) . Z c (“  It  r " 0 T 5 » p,  CO-ODO.M  , .4(1*1  .EPS 
COMMON /G/NST*.  . N,  |TMAy 
INTEGER  ^ . S S 
CEDTP-  G. 0 1 7A*  3 ’91 
C 

C . . .C  ALC'JL  ATE  IN  CLIN  ATTN  ANGt  = S IN  PAOUNS 
C 

the t a( n ■ s.  • 

DO  10  I • 1 . R 

Thc  T a ( l ♦ 1 1 ■ TWFTA(I)  ♦ 10.0  • CCP'T» 

10  CONT I NUF 
C 

C . . ,C  a LC 'JL  AT  F »!I“IJTM  ANGLES  IN  PADIANS 

c 

»i;ni«  io.*cpotp 

DO  20  1-1*17 

phi  ( i*  i ).  20. *C  = DT? *om! ( I i 
20  CONTINUE 
C 

C..  .CALC'JL  ate  ALL  Th  E DTPCCTITN  ESIN'S  QC  SOURCE  SECTORS 

c 

00  A 0 1-1.9 
ZMn-siN(TH=TMin 
DO  ao  J • 1 . 1 *» 

* s ( I.J  )-COS(TurTA( I ) ) • 0 0 * ( 0 h T ( J 1 ) 

t s ( i . j i • cos  t twc  ta  n imsinohium 

AO  CONTINUE 

c 

C ..  .C  ALC'JL  ate  THE  -)T°r*nnN  r0SlNrs  F^c  ThF  N^cmal  vFCTGR  CE  ALL  plan'AP  LEAF 
C INCLINATION  ANGLES  a'S"-!NG  That  THE  a:iwUTh  ANGIE  IS  f OUt  L TO  ZERO  OEGPEFS. 
C 

OH  10  I • 1.9 

XL  E ( I > - -S  IN (Tm:t» ( I ) ) 

YLC ( 1 1 • 0.0 
Z L P C I > - COS(Tucta(  I ) 1 
10  CONTINUE 

c 

C..  .CALCULATE  THE  A»«FVJTe  VALUE  OF  Tnc  n^T  Pan-vJCTS  OF  ALL  S PIJOC  6 -l  f A F 
C ANGLE  p=pyuTAT»  gns.  YhTS  VALUE  IS  EQUAL  TO  THE  COSINE  FACTPF  DESIRED. 

C 

00  SO  L 1-1.9 
00  SO  sb-  i ,o 
DC  SO  SS-  1.1° 

dot-  (UMLn*YS(sn.sf)*vLE(in*YS(sp.ss)*;iF(Lii*:s<st')) 

COSTA(Ll.Se.SS)-  A«s  (PIT) 

50  CONTINUF 
o F TUPN 
END 


SURPOUTIN'C  r.onwoFf 
C 
C 

C . . . SURP  DIIT  I NE  GPONDgy  CALCULATES  THE  Thf  o R a L G90IJN0  FXITANCE  GIVEN  THE  TPUE 
C GROUND  SUPFACF  tfmoc  9 a Tl'»c  . 

C 

C 

Com-pn  n / vw.pm,*l • n c 3 1 

COMMON  /RZ  PGAP(l.9|.  PHI  T(  1,0)  . frlAB?(  3,  q)  , pm  T T 2 ( > . Q ) , $ TF F 
COMMPN fZ /C ns t A ( q, 9 , 1 0 ) , = - T * « V ( A ) , a <» S pa n ( l ) , t««Y,  fGPO*  S E C T A a ( 9 ) 
Common  n /C  QWT  M.  S , o i , c < 1 , S , 9 » , S')MT  ( 3 * Q 9 * <ELV.  GT,  N'JSIK.  I T I ME 
COMMQN/F  / AT.ThF  T A ( 9 ) ,OhT  (l  » ),YLF  (9).  Y L F ( 9 ) , * L F ( 9 ) , Y S ( 9 , 1 «*  ) 

COMMPN/F /YS  (9, 1«  ) . 7 S < ’ 1 . :fdtj,  B,  f p EQn ( 9. 1 ) , WAU5J.EPS 
COMMON  / G/  NSTG , N,  IT-4* 

C 

EGRO*  Ef*ISSV(  A 1#STCF*  (CT*R  ) A 

f. 

«c  T'lPN 
END 


rHIS  p»0£  IS  B&si  au«.m  motuvutt. 
no M OOr*  WBUUSrt®  TO  DDC  
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SUBROUTINE 

C 

c 

C . . .S'JBRQJTI  NE  SKVCif  CALCULATES  ThF  Tmts-aL  F»ITaNCE  FPOM  THE  S<Y. 

C 

C 

Cn****ns  / a / u'» , *>.*, »t  #T)  I **> 

CH-^-nS  /a/  PMT(Ur.).I'r,»‘>?P.<5»,  PhIT?(3,9),STFF 

CQ**rN /C/C05TA  |Q.  9.  1 «*  ) , = - I SSV  (4  ) , abspR®  ( E GpD  » SECT  ap ( 9 ) 

:n>"ipN/'’/rn"TP^.T)»rn,,i,T|,  cn«Tf3»9t,  «flv,  GT  # NUSIH.  I time 
cn',-n/f/u.NPTim,5MinQi,if!.Fm,  vifo).  HFm#  ys(9,18) 
CrOTP.  P,  FRFCK9.31  . '«iA(l51»EPS 

CO^CN/G/NSTS.n,  I t **  A J 

F • 1.0*C.2C>1*'  VP(-?,77C»<,*(AT)**?| 
f SKY -STE F • ( a T* n| • • A • F 
C 

R f TUPS 
END 


SU3PCJT1SS  •*  A 


..SUBROUTINE  CANS?**  CALCULATES  Tmt  CA*>n»v  'JfT'F  Tp  Y C V F c F ! C I E N TS  • 

• • T*<  F S j*-p  G'JT  I NF  CHr.cOM  CALLS  SU*®nHTlNF  So  VC  "PH  WHICH  IS  A *COIFIFD 
PORTION  OF  Thf  SovC  “0 0e  L THAT  CALCULATES  THE  CANOPY  OEP*jETPY 
PARAMETERS. 

CO--ON/GE  0/  OMITU^.P).PLAm.n,SLM(3.n>AYLFAL19,3)»AYLFA(19»3> 

CO'4-  0*1  / A / wv»om,oi*?MM 

CO^PN  f®/  **A0(*»9|»  °HT.  T < ■»  ,9  > , P-.  AO  2 l 3 . 9 > • ®MTT?I  "**91  »STFF 
CHM-PN /C/CH«  T A ( O,  9,  1 a ) , C - I ssv  ( a 1 . a PPHOP  ( •»)  , ESKY,  EGRP,  SECTAP(9) 
co^^rs /o/:ruT( i , a , o ) ,c (y ,« .9) , ?i'«tm.q)v  *-lv.  gt.  nusi*,  hue 

C P*MCN  fl  t\  T . TwC  T A ( 'i  ) » P-^  * ( IP  ) , * L = < 9 ) , Y L c < *3  I • Z L F < 9 J , ‘*>(9.1°) 
CCi-r'I/f/YSn.^h^n),  CCOTO,  a,  F fl  E CO  ( 9 » ? ) , U A ( 1 S ) » t PS 
CC4“CN/G/NSTC,n,  TT**ak 

call  S p vc-oo 

00  10  1*1.3 
00  1C  H.l ,9 

..transfer  IOcnti-al  aopa  vc  phIT  ano  phITI.  pwit  contains  Thf 
opna  ^ 1L  I Tt  OF  w]T  C T = f F I C I f N T S Fpp  EACH  VIEW  ANCLE  ANO  LAYER 
P=  D-.JTA  t pw 

PH  I T(  I , yi  -P^T  Ml  r • “ » 

..CALCULATE  The  °po3  ABILITY  of  GA»  <p-,ap>  for  Ml  ® c« *UT AT l CNS • 

P G A P ( I »“>• 1.-PHJT  f T . *) 

. .CalC'JL  aTF  The  sb-’*5  ILITy  pF  'A9  A vp  hIT  FOR  t^f  HALF  LAvEES(PGAP2#pm[T2» 
FOR  ALL  PFBHJTAtJHNS. 

P5AP?(I.M)-  S?OT(PSAPfI.H>) 

®HIT?(  l #M)  -l.-P^APP  (I.  H) 

10  CONTINUE 


..OBTAIN  The  f cc  oijcsry  of  P^CUPcncf  (ForoO)  Cc  ELEMENTS  IN  EACH  OF  ThF 
..NINE  INCLINATION  ISTCPVALS  FOR  EACH  L A Y F R , 

00  IS  J • 1 » 3 
A 0 0 • 0 • 0 
00  20  N-1.9 

FP=0D(N.JI*  avlca(’*n,j» 

APO.AOO  ♦ FOEOO(N.J) 

20  CONTINUE  . v 

TP  5 5 <.1.9  cV'A-'' 

00(  * > J I »c»f  Ollt.  Jl  UV)  ,\'i  1 

’5  cririMUF  ,4)0 

15  COST  IN  Uf  , ,V  * fO^ 

-ETU.N  V’J4' 


>-^r 
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1 

c 

SUaeO'JTINF  s»vcnpo 

c 

r . . 

• $ t • 0 P Oil  T I N F sev'^CO  T * a “THFirn  vfpfION  Of  ^ °n»TIOkl  r f The  S»VC 

5 

c 

c 

MnnEL  WHICH  CALCUALTFF  the  GFO“FTPir  papa“FTFPS  "F  a canopy. 

c 

co'***nN/r,Eo/  HiTin.oi.ctUf'.n.M  mi  m i.su^h  ]«>,  ? 
•MANSION  NAM*UFn,^i,rt^(A.i,n).TwFTA(101 

) . AYl F A(  19,  3 ) 

10 

c 

OIM^NSirN  Pm;t(  ’M,IOI  .',Tjp),n!)«no,,*t<,q|,rLFM  191 

0 1 MF  NS  tCN  TlfAIlQI,  ( l 7 J . nm.  np(9) 

o=AL  INCLF 

15 

c 

c . . 

• .3E  NEPAL  S Ihih  * T I "N  :hn<;tjaints 

S c VC 

c 

CFrir’»  1 . J 70706  Y ? 

cezpi*  fc.2«n«*3a 

C E 1 » I ■ 3 . 1 A 1 5 9 ? S 5 

?0 

Ct  0 T*«  . 0 1 7A*  » ?<?3 

C F 3 T 0 • 57.?<35‘r7‘3‘5 

Cfwtc.  . CCG  ? RO  a P fl  2 1 
naan  ns  • <? 

N-  AT  • 1 

2 5 

r 

c . . 

N L A Y . 3 

9 A N 0 W • 3 C / N 9 A N 0 S 

S e VC 

. .PAB a^etfp  ini  7! £i I Nation  and  convepsion 

s®vc 

30 

NSTje-N9AND<  *-T 

S 0 VC 

r 

PANDw • flANQy*CFnT0 

SPVC 

L 

r 

..  C PE  = C i:  IF  NTS  F n o nicFMSF  aurAHCV  V-CTTPS 

SPVC 

35 

V. 

AL  pH  A?  ■»  0. 

S*VC 

SINA--0. 

s P VC 

00  ? I-l.N^ANHS 

s °vc 

SI N Al • S INA? 

S p VC 

ALpHA?«  *LPm»,*°AN')W 

SPVC 

<■0 

S I N A ? • s I N ( t 'w  0 u A ? I 

SPVC 

*K  fl  ) "S  J A A?*S  TNA'*-.^INAl*STNAl 

save 

? CONTINUE 

c 

c . . 

..SOURCE  3IP?CMPM  INCLINATION  ANGLES 

F-tVC 

<.5 

c 

T” T a l .0. 

SRVC 

thet  a i : M^'nw/?, ) *9ANnv 

SPVC 

on  * l"l»NHANHC 

S F VC 

THE  T a ( I*1J«thctmi  )*m‘nu 

SPVC 

50 

r 

C. . 

3 CONTINUE 

..CANOPY  Rffl-fTOV.  P A r M f A N H P Y [ »vro  IS  CO“3nPCD  OF  0N£ 

OPTICAL 

cp  vc 

C . . 

..“AT-PIAL  - “AY  *C  SPEC  I FTC*.  AND  UN  J CUE  OFO“ctsICal 

pp  n°EF TIES. 

FB  VC 

C.  . 

..CAN03Y  Gf  f'*-*1  T = jr  PA5A“fY  = PS  CONST  FT  Of  ( 1 ) L c A F A Nr,  It  fOfCUENCY 

ssvr 

55 

c . . 

. .0  IS  TPI9UTICN  CIIN-TIPN  nCNOTCO  PY  X l F A AND  YLFA  (?)LfAf 

A p F * I NOE  X 

s°vc 

c . . 

NOTED  PY  fiat  and  ( 3 ) c A N o 0 Y ^ F N S I T Y 0 F N p T f o «v  s l M . 

Y l E A ( HFC ) 

A!  \JC 

c. . 

..AND  YLFA  KJST  PE  <PFCICTCD  AT  AN  TOn  NU**opa  (NANO)  CF 

EVENLY  SPACED 

SPVC 

c .. 

r 

..POINTS.  FLAI  IS  now-wfga Tt VE  and  SlAl  RANGES  BETWEEN 

0 AND  1. 

SPVC 

60 

v» 

OF  L F • 10. * C E nTo 

SPVC 

00  350  U-l.NLAY 

SPVC 

jCHiS  FASt  IS  BiSl  ftUALin  PKACUUW* 
raoM  ooi-  'i  ?u&uj..'>fUiL  ro  — 
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65 


70 


75 


eo 


85 


90 


95 


100 


105 


1 10 


115 


120 


125 


no 


NANG*  1 9 

C 

C . . . ASS  I GN  THE  NUm3c  R nc  *«ATE<IA15  IN  ANY  GIVEN  l A Y E 0 
C 

I M A T ■ 1 

-Tp( u ) • 1 “ a T Save 

I - A T 1 • I *A  T 
00  551  J ■ 1 1 1 * A T l 

n a t - j 

p P <.  1 -M-1,NANG 

*L  f & m-  ) ■ a <LC  A ( 1 1 ) 

Y L E A ( - - ) ■ A Y l P A (Mu*  It  ) 

5 1 CONTINUE 

c 


C ....  INTEGRATE  an-)  uramEiye  r^-  LrAc  ANGLE  ERCO'JFNCY  DISTRIBUTION  5aVC 

C . . . . P Un:  T n*-'  USING  5T-E >S0n<  5"Lr..Tm<  IS  DENOTED  5Y  E.  5avC 

C....“-l  -DUALLY  ISTC9VAI.S  OF  f ARE  Thcn  OET=R-.INED  AND  DENOTED  SRVC 

C....5Y  TLA  <-  e C T N T G ) , THE  TA^LF  FLA  IS  USED  FOR  RANOfMLY  SELECTING  SRVC 
C....LCAF  INCLINATION  ANGL  F * • $ = VC 

r 

DO  505  !«1,van G SRVC 

3C5  Xlfill  )»XLF.*.{  M*'.C}T3  S = VC 

-«(  (NANG-1  )/’  Ml  «5VC 

NANSLE ( II • Im'  T ) «m  S "VC 

CALL  Tal?(’4tUcA>'ri:i,n’<lC|  S“VC 

DC  310  I A NG» 1 • M S»VC 

310  EL A(  IL » I* AT. I AMr» j,pM( JANG ) S 9 VC 

C 

C . . . . NQRMA  L I zr  THc  T NO  I]T  L c A p EREOJ^N^Y  DISTRIBUTION  FUNCTION  TO  C s- T A I N 53VC 
C....A  DENSITY  FUNCTION  f WHICH  IS  SPECIFIED  AT  ►*  “TINTS.  SRVC 

C 

FT  or. 3,  grvC 

DO  311  I-l.NANG  S^VC 

311  ETOT -c TOT  * YL c A ( T ) S“VC 

DO  312  1-1.9  S=VC 

312  F c t) - ( YL l A ( ’* T 1 *YL c a ( ? n If  TOT  S?VC 

DC  315  l « 1 , n * ng  S«VC 

315  <LFA{I)»XLrA(I)B-COTO  53  VC 

*-M*l  S 9 VC 

...CALC  JLATF  Tt-E  -s\m  t>  o 0J«  C T T ’“N  frt>)  IN  The  OI=*tCTION  nr  The  SOURCE  5 o vc 
C . . . ( THE  T a ) OF  ONC  IJNIT  ~ A c A3=ft  uTTm  INCLINATION  INClF.  ToE  LEAVES  SRVC 
C....AT  thj$  AKGLr  A»E  s'S'l-ED  T H «c  A 7 1 J T-»  AL  L V ISOTRPPIC.  SRVC 

C 

DO  ’33  r ANTI f *2r V$nM©  SRVC 

INCLF*  -5.  *:n  T3  SFVC 

DO  3 2D  1-1.3  53V c 

INCl F • l* CLE -oclc  SRVC 

320  CALL  COP C INCL F. Tmct A ( T ANGLE ) , re ( n ,cpp  102) 

C 

C... .CALCULATE  T Hc  -FAN  PROJECTION  (?““)  IN  TmE  DIRECTION  pp  THE  SOURCE  SRVC 
C • • • • ( The  T A ) Cr  ON?  "NIT  lfaf  AREA  AVERAGED  OVER  ThE  CANOPY  LEAF  ANGLE  5 R VC 
C...  .DENSITY  FUNCTION  c.  SDVC 

C 

call  cnp-ir.oo.oPMf tangle n save 

c 

C. . . .C  ALCULATF  T*c  op  0*5  ABILITY  OF  A NTT  (PHIT)  coo  a light  RAY  WITH  SRVC 

C... .SOURCE  DIRECTION  thRTA.  SRVC 


c 

CALL  POENS(TL.TMAT,lANC.LE.00>MUk,OLc)tTHET\.NANGLF>ELA,SlAl»FLAl> 
• PHIT) 

330  CONTINUE 
35  ] C ONI  j NiJE 
350  CONTINUE 

J -N- at 

DO  1-1,3 

or  ’2*3  i fO 

PH  IT  1 I I 0 - | -OMT  T ( | , J , M*|  | 

2?B  CONTINUE 
RETURN 
ENO 


this 

tfKOS* 


i ^ 


io  bdc 
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1 


5 


10 


15 


20 


SU^BC’JTINt  CO®ML°MA»°c  TA,0®,rcPJn>) 
r 
C 

C....TvMS  PHnG®A.M  ;hC"UtF^  t*<'  “f  an  pea jc' tj  rv  of  \ i»njt  leaf  abFa  IN 
C....THE  DIRECTION  fF  Thf  5 ni'orc4  TMr  i c a P js  I*i"LIw£0  AT  an  ANGLE 
C...  .ALPHA  A*:n  is  aSSM«c}  IT  A,MMUTMALIY  ISPT3TMC. 

C 

c 

c 

op »c ns  ( ai  p h a > *m«i^tai 

IP  (Al  PHA.LF  .«=  T.\)  DST'.ICN 

C 

C • • • • The  T AO  IS  Tn-  If  AC  4>Tu!!n  ANGLF  AT  U»4l  CM  f13  ^ECC'Tf*  NEGATIVE  AND 
C....IS  Is1  THE  f l p f t Q".\r>HANT.  Thf  FUNCTION  n P I G ^y^-ETBIC  AN’O  HENCE 
C....IS  AVfPAGFO  TVPP  Lc AP  A7I-UTM  ANGLES  CF  0 TO  PI  CaCIANS. 

C 

THEUO- ACnSfM^MrTSJ/TANfAlPWM) 

TANTO* T AN ( Thc  T aO  J 

pp«PP*(l.»(Tl*JTO«THPTAC)/rFPir?) 

P F T U ® N 
PN0 


r op 


c op 
cop 

cop 


cop 

COP 

COP 

COP 

COP 

COP 

CCP 

COP 

COP 

COP 


1 


5 


10 


15 


SU3®  O'  I TINE  roarer. , pa.np'ii 
C 
C 

(••••This  P*nG®A“  GALOMIMPS  th®  m=»n  pooJECTTCN  Or  a l 1 n I T LEAF  ape  a in 
(••••The  DIRECTION  nr  Twf  Sp’orr  (Pd,4>  FT®  TwF  SITUATED  CANODY.  THE 
C... .LEAVES  r F Thc  CANH3Y  a»f  aSSUmfh  Tn  *F  aZINMTmaLLY  ISPT®CPIC.  ThF 
C ....  OP  cijn:  TICN  USFO  TN  Th®  : a t C »‘L  a T I r n MAS  3 F E n PP  = ’.'ICUSLY  OtTcffMINEO 
C....FQR  A GIVEN  Sn-IPCE  OTPFCTION  F no  c A F INCLINATION  ANGLES  OF 
C....5.  15,  ...»  «5  OEG^EFS.  G IS  THE  LEAF  INCLINATION  ANGLE  DENSITY 
C... .FUNCTION. 

C 

C 

ci  v:*  s rc\  cprn.cip  ) 

OP*  -o. 

00  1 I ■ 1 » 3 

l op^,ppH*cp<  n *r,(  n 

PFTUPN 

ENO 


CP®" 


C r pm 
r<3pM 
C'lPN 

CC  P‘1 
C fPK 
C.QPM 
COP" 
CCPH 

c r®  ^ 
r o®  h 
r 3®" 
C D0M 
COp" 
C0®“ 


1 


5 


10 


15 


20 


25 


30 
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